Interaction corrections at intermediate temperatures: Longitudinal conductivity and 

kinetic equation 
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It is well known that electron-electron interaction in two dimensional disordered systems leads to 
logarithmically divergent Altshuler-Aronov corrections to conductivity at low temperatures {Tt <C 1; 
T is the elastic mean-free time) . This paper is devoted to the fate of such corrections at intermediate 
temperatures Tt > 1. We show that in this (ballistic) regime the temperature dependence of conduc- 
tivity is still governed by the same physical processes as the Altshuler-Aronov corrections - electron 
scattering by Friedel oscillations. However, in this regime the correction is linear in temperature; the 
value and even the sign of the slope depends on the strength of electron-electron interaction. (This 
sign change may be relevant for the "metal-insulator" transition observed recently.) We show that 
the slope is directly related to the renormalization of the spin susceptibility and grows as the system 
approaches the ferromagnetic Stoner instability. Also, we obtain the temperature dependence of the 
conductivity in the cross-over region between the diffusive and ballistic regimes. Finally, we derive 
the quantum kinetic equation, which describes electron transport for arbitrary value of Tt. 

PACS numbers: 72.10.-d, 71.30.-|-h, 71.10.Ay 
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I. INTRODUCTION 

Temperature dependent corrections to conductivity 
due to electroiigelectron interactions Jias been a subject 
of theoreticalEl'El and experimentafrEj studies for more 
than two decades. Recently the interest iiL the matter was 
renewed with appearance of new dataOa showing a sign 
change in the temperature dependence of conductivity in 
two dimensions (2D). Theoretical discussionala that fol- 
lowed emphasized the question of whether that data indi- 
cated a non-Fermi liquid behavior. However, the experi- 
ments were performed in a regime where the temperature 
T was of the same order of magnitude as the inverse scat- 
tering time (obtained from the Drude conductivity), 



while pre-existing calculations were, focused on the two 
limiting cases: Xhc diffusive regimeol 
ballistic regimeDu Tt ^ 1. 



Tt < 1, and the 



limiting cases: Xhc diffusive regime 
allistic regimeDu Tt ^ 1. ^ 
In the diffusive limit one findsiji^ for the logarithmi- 
cally divergent correction to the diagonal conductivity 
Sa: 



da = 



27r2;i 



In 



h 

TV 
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(1.1) 



where Fq is the interaction constant in the triplet chan- 
nel which depends on the interaction strength. It is clear, 
that the sign of this logarithmically divergent correction 
may be positive (metallic) -or negative (insulating) , de- 
pending on thg value of F(5^E3. 

The rcsultP'Ll for the ballistic region frequently cited in 
literature reads 



Scr = 



e 



Tt 



firs), 



(1.2) 



where /(r^) is a positive function of the gas parameter 
of the system, r^. In a sharp contrast to Eq. (1.1), equa- 



tion (1.2) predicts always metallic sign of the interaction 
correction. 

The absence of a rigorous calculation at intermediate 
temperatures, Tt/Ti ~ 1, may have contributed to the 
notion that those two limits are governed by different 
physical processes. In this pap er w e prove that notion 
erroneous: the results (1.1) and (1.2) are due to the same 
physical process, namely elastic scattering of electrons by 
the self-consistent potential created by all the other elec- 
trons. Therefore, these two different expressions are in 
fact the two limits of a single interaction correction. We 
calculate the correction within assumptions of the Fermi 
liquid theory (other limitations of our approach we dis- 
cuss below) and present the cross-over function between 
the diffusive and ballistic limits. 

Moreover, w e sh ow that the existing theory for the 
ballistic limit (1.2) is incomplete. First, the results of 
Ref. account only for Hartree-like interaction terms 
missing the exchange or Fock terms. Second, this theory 
essentially employs a perturbative expansion in terms of 
the interaction strength, which breaks down for stronger 
coupling. Both issues lead to the change in the theoreti- 
cal prediction even on a qualitative level. 

The consequence of the first point is that the correction 
to conductivity (1.2) is always negative unlike the correc- 



tion in the diffusive limit that changes sign depending on 
the value of . This sign change is due to competition 
between the universal (and positive) Fock correction and 
the coupling-specific (and negative) Hartree contribution. 
If the Fock potential (or, to be more precise, singlet chan- 
nel) is properly taken into account, then the sign of the 
correction in the ballistic limit is also not universal (being 
positive for weak interaction in contrast to Refs. ^-0). 

As follows from the second point, for the stronger in- 
teraction the Hartree correction should be modified to 
include higher order processes. For this case we show 
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(see Section III ) that in fact it should be replaced by the 
triplet channel correction, which is characterized by the 
Fermi liquid constant . This constant measures the 
strength of the spin-exchange interaction. If < 0, the 
spin-exchange interaction tends to align electron spins 
and (if it is stroug, enough) leads to the ferromagnetic 
Stoner instabilityc^l. Even though this constant is un- 
known, it can be found experimentally by means of in- 
dependent measurement of the spin susceptibility of the 
system. As a function of temperature the interaction 
correction to conductivity is almost always monotonous, 
except for a narrow region of parameters (where it is so 
small that it can hardly be observed). 

The remainder of the paper is organized as follows: the 
following section is devoted to qualitative discussion of 
the physics involved. In the same section we summarize 
our results. Then we present two alternative approaches 
to the microscopic calculation. Jn Section |l| we use the 
traditional perturbation theoryEI to derive the results pre- 
sented in Section while in Section IV the same results 
are obtained using the quantum kinetic equation that we 
derive. The advantage of the kinetic equation approach 
is that it can be readily used to discuss the temperature 
behavior of quantities other than conductivity. These re- 
sults are advertised in Conclusions and will be published 
elsewhereEJ. 



Consider a single impurity localized at some point, 
taken as the origin. The impurity potential U{r) induces 
a modulation of electron density close to the impurity. 
The oscillating part of the modulation is known as the 
Friedel oscillation, which in 2D can be written as 



Sp{r} = - 



27rr2 



sin(2fci?r) 



(2.1) 



Here r denotes the distance to the impurity and its poten- 
tial is treated in the Born approximation ^ — J U{r)dr. 
In 2D the free electron DoS is given hy v — mjixh^ and 
TO is the electron mass, kp 1% the Fermi momentum. 

Taking into account electron-electron interaction 
V^(f\ — r2) one finds additio nal s cattering potential due 
to the Friedel oscillation Eq. ( |2.l| ) . This potential can be 
presented as ja^sum of the direct (Hartree) and exchange 
(Fock) termsllZl 

bYif^^r^) = VH(fi)(5(fi - fa) - VF{rx.r2)\ (2.2a) 



(2.2b) 



(2.2c) 



II. QUALITATIVE DISCUSSION AND RESULTS 

In this section we describe the scattering processes con- 
tributing to the temperature dependence of conductivity. 
We show that unlike the standard Fermi liquid correc- 
tions, the leading correction to conductivity is accumu- 
lated at large distances, of the order VFl'min{T, yjT/r). 
In the ballistic limit such correction is linear in temper- 
ature and we derive this result here using a text-book 
quantum mechanical approach. The diffusive limit is 
discussed in detail in Ref. |l]. The resulting correction 
5a InT seems to be rather different from the linear 
one, but we show that both corrections arise due to the 
same physics - coherent scattering by Friedel oscillations. 
Throughout the paper we keep the units such that h — 1, 
except for the final answers. 



A. Scattering by Friedel oscillations 

We start with the simplest case of a weak short-range 
interaction Vo{ri — fa) and show how one can obtain the 
correction to conductivity in the ballistic limit, i.e. due 
to a single scatterer. This discussion is similar to that 
of Ref. where the correction to the one-particle den- 
sity of states (DoS) was discussed, and also of Ref. [l^ , 
which describes the correction to the conductivity in the 
diffusive limit. 



where by p{r) we denote diagonal elements of the one 
electron density matrix n. 



l(^l,f2) =^«'fe(^l)«'fe(r2). 



(2.3) 



The factor 1/2 indicates that only electrons with the 
same spin participate in exchange interaction. As a func- 
tion of distance from the impurity the Hartree-Fock en- 
ergy 6V oscillates similarly to Eq. (2.1). 




FIG. 1. Scattering by the Friedel oscillation. Interfer- 
ence between the two paths A and B contributes mostly 
to backscattering. The Friedel oscillation is created due to 
backscattering on the impurity, path C. 
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The leading correction to conductivity is a result of in- 
terference between two semi-classical paths depicted on 
Fig. |l| If an electron follows path "A" , it scatters off 
the Friedel oscillation created by the impurity and path 
"B" corresponds to scattering by the impurity itself. In- 
terference is most important for scattering angles close 
to TT (or for backscattering) , since the extra phase fac- 
tor accumulated by the electron on path "A" (e*^'^^ with 
R being the length of the extra path interval relative to 
"B" and 2k being the difference between initial and fi- 
nal momenta for that extra path interval) is canceled by 
the phase of the Friedel oscillation e"'^*^^^ so that the 
amplitudes corresponding to the two paths are coherent. 
As a result, the probability of backscattering is greater 
than the classical expectation (taken into account in the 
Drude conductivity). Therefore, taking into account in- 
terference effects leads to a correction to conductivity. 
We note that the interference persists to large distances, 
limited only by temperature i?^ l/\k — kp\ < vp/T. 
Thus there is a possibility for the correction to have a 
non-trivial temperature dependence. The sign of the cor- 
rection depends on the sign of the coupling constant that 
describes electron-electron interaction. 

To put the above argument onto more rigorous footing 
and to find the temperature dependence of the correc- 
tion con sider now a scattering problem in tha potential 
Eq. (2.2). Following the textbook approachlla, we write 
a particle's wave function as a sum of the incoming plane 
wave and the out-coming spherical wave (in 2D it is given 
by a Bessel function, which we replace by its asymptotic 
form) 



^ = e 



ik-f 



(2.4) 



Here f{9) is the scattering amplitude, which we will dis- 
cuss in the Born approximation. For the impurity po- 
tential itself the amplitude f{9) weakly depends on the 
angle. At zero temperature it determines the Drude con- 
ductivity ao, while the leading temperature correction is 
proportional to T^, as is usual for Fermi systems. We now 
show that this is not the case for the potential Eq. (2.2). 
In fact, taking into account Eq. ( ^.2| ) leads to enhanced 
backscattering and thus to the conductivity correction 
that is linear in temperature. 

First, we discuss the Hartree potential Eq. (2.2b). Far 



e can 



from the scatterer the wave function of a particL 
be found in the first order of the perturbation theoiy as 
yj, _ gifc r _|_ |5vl/(r), where the correction is given bylla 



(5*(f) = i / dfiVH{fi)e 



ik-ri 



2tt 



k\f— fi I 



ik\f—fi I 



(2.5) 



Here |r — ri| ~ r — r • ri/r, since we are looking at 
large distances. Substituting the form of the potential 
Eq. (2.2b) and introducing the Fourier transfer of the 
electron-electron interaction Vq we can rewrite Eq. (2.5) 
as 



(5*(r) 



where 



/2tt 



ikr 



^sin(2fcfri)e'«-""\ (2.6) 
kr J rl 



q = k-kf/r, \q\^2ksm{e/2), 

with 6 being the angle of scattering. Comparing to 
Eq. ( |2.4| ) we find the scattering amplitude as a function 
of 9 (it also depends on the electron's energy e = k"^ /2m) 

f{0) = -^^o(g) J ^ sin(2fc;^r)e'^-: (2.7) 

The integral can be evaluated exactly0 and the result is 
given by 



f{0) 



Vo{q) 



\q\ < 2kF; 



277 I arcsin(^) , \q\ > 2kF. 



(2.8) 



Let us examine this expression more clo sely. 
\q\ < 2k, the scattering amplitude Eq, 



(2. 



Since 
for small k 



weakly depends on the angle through the Fourier com- 
ponent of the interaction Vo{q), see background value of 
f{9) on Fig. ||. However, we are dealing with electronic 
excitations close to the Fermi surface, so in fact k is close 
to kp, \k — kp\/kp ^ 1. If k > kp, then the scattering 
amplitude Eq. (^|^) has a non-trivial angular dependence 
around 9 ~ tt shown on Fig. ^j. 

According to Eq. ( |2.8| ) such dependence is only pos- 
sible in the region |g| > 2kp. This translates into the 
condition 16* — 7r| < [2{k — kp) /kpY/"^ , which determines 
the singular dependence of the width of the feature in 
the scattering amplitude on the energy of the scattered 
electron. Finally, using the fact that arcsin(l — x) = 
-K 12 — \/2x, we find that the dependence of the height of 
the feature in the scattering amplitude is also singular: 
8j{d)^\{k-kp)lkp\^l\ 



m 



71 




e 



(£-£ /£ )l/2 
F F 



FIG. 2. Scattering amplitude. The singularity for 
backscattering is due to interference of paths "A" and "B" 
on Fig. 1. 
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The transport scattering rate is determined by the 
scattering cross-section and can be found with the help 
of the amplitude Eq. (2.8), as well as the constant am- 



plitude fo of the scattering by the impurity itself 



^(l-cos0)|/o-t-/(e)r 



(2.9) 



The leading energy dependence of r^^ comes from the 
interference term, which is proportional to f{9). Then 
integration around 9 = tt is dominated by the feature of 
f{9) resulting in a term of order (e — eF)/f-F- It is this 
term that gives rise to the linear temperature dependence 
we are after. Since we are interested in this leading cor- 
rection only, in all other terms we can set k ^ kp and 
write the scattering rate as 



z tp 



n{t-tp)h. (2.10) 



Here r\(x) is the Heaviside step function and t^^ is the 
zero-temperature rate that determines the Drude con- 
ductivity (indeed, the 9 ^ tt feature in f{9) only exists 
for k > kp and at T = there are no electrons with 
k > kp). 

To obtain the scattering t ime w e have to integrate the 
energy-dependent rate Eq. (2.10) with the derivative of 
the Fermi distribution function np{e) 



f d 
r = / deT{e)—np{e). 



Then the second term in Eq. ( 2.10| ) leads to a linear cor- 
rection to the Drude conductivity, small as T/ep. How- 
ever this is not the only contribution to the temperature 
dependence. At finite temperatures we also have to mod- 
ify the Friedel oscillation Eq. (2.1) as follows: 



5p{r) = - 



2TrVp sinh^ 



sin{2kpr). 



Consequently, the scattering amplitude Eq. (2.7) be- 
comes temperature dependent 



fiO) = -"^Voiq) 

ZTT 



^ e^'^'''~ sm{2kpr2)e'^-^' 



S T 
— = ^2i^Vo{2kp) — . 



(2.12) 



The conductivity Eq. ( 2.12 ) is the same correction as 
the one calculated in Ref. |6| see also Eq. ( |l.2D, up to 
a numerical factor. It is also clear that Eq. ( |2.12 ) is 
not the full story. We ha ve f orgotten about the Fock 
part of the potential Eq. ( ^^ ! Substituting Eq. ( |2^ ) 
into Eq. (2.3|), we find the perturbation of the d ensit y 
matrix [which appears in the Fock potential Eq. (2.2c)] 
(5n(ri,r2) ~ Sp[{ri -\- r2)/2]. Then the argument can be 
repeated. The only difference is that the leading tem- 
perature correction comes from the Fourier component 
at g = 0, rather than q = 2kp. What is most important, 
the Fock potential enters with the opposite sign. There- 
fore the expression for the conductivity Eq. (2.12) has to 
be corrected 



a = ao 



{2Vo{2kp)-Vo{Q))^ 



(2.13) 



The sign of the correction is thus not universal and 
depends on the details of electron-electron scattering. 
If the weak interaction is reasonably long-ranged, then 
Vb(0) > Vii{2kp), so that the correction in Eq. ( |2.13| ) 
has the sign opposite to that in Ref. 0. 




/^A . f^kp 

-T^Vb (7)arcsm 

27r \ P 



(2.11) 



P=y(^) +(9 + 2M' + y(^) +{q-2kpY. 

Neglecting the small temperature dependent term in the 
denominator in Eq. (2.11) brings us back to Eq. (^.8]). 
Keeping this term leads to the same feature in f{9) as 
the one on Fig. ^, only now its width and magnitude are 
proportional to -v/T. The resulting correction to the con- 
ductivity is therefore similar to the one discussed above. 
Up to a numerical coefficient. 



FIG. 3. Scattering process with two impurities and the 
Friedel oscillation. Scattering to all angles is affected by in- 
terference. The relevant Friedel oscillation is created by the 
self-intersecting path C. 

So far we have considered the effect of a single impu- 
rity. The extension of the above arguments to the case 
of many impurities is straightforward. In particular, one 
can consider a scattering process, which involves two im- 
purities and the Friedel oscillation shown on Fig. |[ It is 
clear that this process contributes to the scattering am- 
plitude at any angle, and not just for backscattering as 
the single impurity process on Fig. |l| (which is typical for 
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the diffusive motion of electrons). Such processes were 
discussed in detailfalthough using a sUghtly different 
language) in Ref. |6[ Scattering by Friedel oscillations 
created by multiple impurities results in a conductivity 
correction (1.1) that is logarithmic in temperature and is 
typical for 2D diffusive systemsEl. 

Comparing the scattering processes on Figs. and 
^, one can clearly see that conductivity corrections, 
which arise from these processes are governed by the 
same physics: coherent scattering by the Friedel oscilla- 
tion, which means that the ballistic and diffusive regions 
should be analyzed on the same footing. In the next sub- 
section we present the results of such anlysis, postponing 
the actual calculations until Sees. Ill and IV. 



w{x > 1) 



Sttx 



ln(2x) - - (In X - 1) (6 In 2 - 1) 



In the opposite limit Tt ^ 1 it approaches a constant 
value (C ~ 0.577 ... is the Euler's constant and Ci^) is a 
derivative of the Riemann zeta function) , 



w{x < 1) w 1 + — — In x - In 2 



2t:x 
9" 



C + -+6C'(2) 



so that the linear correction does not completely vanish 
in the diffusive limit, but competes with the logarithmic 
term and in semiconductor structures with low Fermi en- 
ergy it might be important except for the lowest temper- 
atures. The full function z«(x) is plotted on Fig ^. 



B. Results 



Let us first consider the case of a weak, short range in- 
teraction potential. Then the interaction can be treated 
in the lowest order of perturbation theory, so that the re- 
sulting correction is proportional to the interaction con- 
stant: 



Saw — 



Tt 



71- 



1 



-w(Tt) 



72 Ep 
An ^ T 



(2.14) 



Similarly to Eq. (2.12), it can be written as a sum of 
Hartree and Fock contributions (similar expression for 
7i in one-dimensional systems was obtained in Ref. |2C|): 



71 = VoiO) - 2Voi2kF) 



12 V^{Q)-2{Vo{k)) 



FS 



(2.15) 



where (. . .)fs stands for the average over the Fermi sur- 
face. Here we kept the notation for the electron-electron 
interaction adopted in the previous Section. Then the 
Hartree correction is proportional to the Fourier compo- 
nent of Vb((7) at g = 2fci?, while for the Fock correction 
(7 = 0. The two corrections have different sign as we 
discussed above. The extra factor of 2 in the Hartree 
correction is due to electron spin degeneracy. 

Note, that Eq.( ^.14 ) is defined only up to a 
temperature-independent constant which is determined 
by the ultraviolet contribution. We have chosen the ar- 
gument of the logarithm to be Ep/T instead of the usual 
1/Tr to emphasize that contrary to the naive expecta- 
tions the logarithmic term persists up to temperatures 
much larger than 1/r, see also Ref. ^ 

The different expressions for the Hartree terms in 71 
and 72 are related to the fact that the single impurity 
scattering, see Fig. |l|, and multiple impurity case, see 
Fig. ||, allow for different possible scattering angles. 

The dimensionless function w{Tt) describes the 
crossover between ballistic and diffusive regimes. In the 
ballistic limit Tt 1 it vanishes as 



10 12 



FIG. 4. Dimensionless function w(x), which is defined so 
that w(0) = 1. 

If the Coulomb interaction is considered, then the low- 
est order in interaction is not sufhcicnt since for the long 
range interaction i'VQ{q « 0) ^ 1. Although the in- 
teraction itself is still independent of the electron spin, 
summation of the perturbation theory depends on the 
spin state of the two electrons involved. In the first order 
correction Eq. (2.14) all spin channels gave identical con- 
tributions. The total number of channels is 4 and they 
can be classified by the total spin of the two electrons: 
one state with the total spin zero ("singlet") and three 
states with the total spin 1 ( "triplet" ; the three states 
differ by the value of the z-component of the total spin). 
For long range interaction the perturbation theory for the 
Hartree correction singlet and tripl et c hannels is differ- 
ent. It is knownld'El (see also Section HI), that the singlet 
channel contribution should be combined with the Fock 
correction as a renormalization of the coupling constant. 
However, the final result is universal due to dynamical 
screening: the singlet channel modification of the cou- 
pling does not affect the result. What remains of the 
Hartree term is the triplet channel contribution, which 
now depends on the corresponding Fermi liquid constant 
Fq. Thus, the total correction to the conductivity can be 
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written as a sum of the "charge" (which combines Fock 
and singlet part of Hartree) and triplet contributions 

o = ajj + 6(tt + 5<Jc, (2.16a) 

where the charge channel correction is given by 



5a. 



e2 Tt 



c 



1 



-J{Tt) 



ttTi h 

and the triplet channel correction is 



— In . 

2TT^h T ■ 



(2.16b) 



5ai 



e2 Tt 

-3(1 



^ ln(l + F-)^ ' 



1 



(2.16c) 



here the factor of three in the triplet channel correc- 
tion Eq. ( ^.16c ) is due to the fact that all three com- 
ponents of the triplet state contribute equally. We re- 
iterate that the corrections Eqs.( ^.l^ ) are defined only 
up to a temperature independent (however not neces- 
sarily Fermi-liquid c onsta nt independent) term, see also 
discussion after Eq. (2.14). 



in the first line of Eq. ( 2.16c| ) and 



1 



:ln 



%/2+ v/2~~ 



1 r. 



arctan \\ \r1 — 1, > 2 



in the second line. For Ts ^ 1 our replacament is well 
justified even wit hin w eak coupling scheme. 

Similar to Eq. (2.14) the dimcnsionless functions ][x) 
and t(x\ Fq) describe the cross-over between ballistic and 
diffusive limits. They are plotted on Fi gs. |^ and ^ and 
full expressions are given by Eqs. ( 3.35 ). The universal 
function f(x) has the following limits 

/(x> 1) w 5-f2(lna;- l)ln2- -ln(2a;)'); (2.17a) 

3ttx \ 2 / 



/(x < 1) w 1 - 7ia;+ -xlnx; (2.17b) 



71 =-^ + ^(c + ^ln2)« 0.7216; 



J I I I I I I I I I I I I I L 




2 4 6 8 10 12 14 5 



FIG. 5. Dimensionless function f{x), defined so that 
/(O) = 1 

We should warn the reader here, that we describe the 
interaction in the triplet channel by one coupling con- 
stant Fq . For the weak coupling limit, it corresponds to 
the approximation Vb(2fci?) ~ (Vb(fc))FS- This approxi- 
mation overestimates the triplet channel contribution to 
the ballistic case for Ts = V^e^ /{nhvp) <C 1- However, in 
this limit contribution itself is much smaller than the sin- 
glet one. For better estimates in this regime one should 
use 



FIG. 6. Dimensionless function t{x,Fo) defined so that 

i(o,Fcn = i. 

The function t{x;FQ) depends on the coupling con- 
stant and therefore its asymptotic form also depends on 
Fq. For very small x ^ 1 + Fq the asymptotic form is 



tix « 1 + Fo") « 1 - 722; + I-xlnx(^3+ 



72 



C(2) 



1 



7r73 



9(1 + Fq- 



(2.18) 



1 



V2 



+ 



18 



C 3 



1 

7+Ff 



ln2 1 
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73 = 1 - 



5Fo - 3 
1 + 



ln(l + F° 



Notice that at Tr —i- 0, Eqs. (2.16) reproduce the known 
resuh (1.1). Let us point out that for numerical reasons 
contributions of scahng functions w, /, t change the re- 
sult only by few percents and they can be neglected for 
all the practical purposes. 

No tice t hat while the charge channel correction 
Eq. (|2.16bD is universal, the triplet channel correction 
Eq. (|2.16cj ) is proportional to Fq , which might be nega- 
tive. That leads to the concl usion , that the overall sign 
of the total correction Eqs. (2.16) depends on value of 



F^: it can be either positive or negative, see Fig. ^. 




-50.0 



0.0 1.0 2.0 3.0 4.0 5.0 

Tr/h 

FIG. 7. Total interaction correction to conductivity. The 
divergence! at Tt/Tl — > is due to the usual logarithmic 
correctiontJ. Curve Fq = corresponds to the universal be- 
havior of completely spin polarized electron gas. The cor- 
recti on is defined up to a temperat ure in dependent part, see 
Eq. (3.32) and discussion after Eq. (U.14 



Combining together all of the above results we plot the 
total correction to the conductivity on Fig. ^ for differ- 
ent values of Fq . The divergence at low. temperature is 
due to the usual logarithmic correctioin. Although the 
exact value of Fq can not be calculated theoretically (in 
particular, its relation to the conventional measure of the 
interaction strength, Ts, is unknown for > 1), in prin- 
ciple it can be found from a measurement of the Pauli 
spin susceptibility 



X = 



1 



F-' 



(2.19) 



where the density of states v should be obtained from a 
measurement of the specific heat (at ^ T <C Ep). 
The constant Fq is the only parameter in our theory 
which describes all the data, including the Hall coefhcient 



and the magneto-resistance in the parallel field. The the- 
ory for interaction corrections in the magnetic field will 
be addressed in the forthcoming papeia. 

The correction in Fig. |^ is almost always monotonous, 
except for a narrow region —0.45 < Fq < —0.25. A typi- 
cal curve in this region is shown in Fig. ^. Note, however, 
that the overall magnitude of the correction in the range 
of Tr in Fig. ^is so small that it can hardly be observed. 



da X 




0.2 0.3 0.4 0.5 

Tr/n 

FIG. 8. The non-monotonous correction to conductivity. 
Note the difference in the overall scale relative to the previ- 
ous figure. 

When the interaction becomes so strong that the sys- 
tem approaches the Stoner instability, Fq ceases to be a 
constant a nd be comes momentum-dependent. Thus the 
resuh Eq. (|2T|) is no longer valid. Although the sim- 
ple condition Sax < suggests that this happens at 
T a 
tion 



(1 + Fq)Ep, the more detailed analysis (see Sec- 
III E) shows that it happens much earlier. In fact. 



the approximation of the constant Fq is valid in the pa- 
rameter region defined by the inequality 



T 



<(1 + K)' 



(2.20) 



see Section III E for the origin of this inequality. We were 
not able to make a reliable calculation of Sa{T) at higher 
temperatures. 



III. PERTURBATION THEORY 



In this section we show how the announced results 



Eq. ( 2.16 ) can be obtained with the help of the tradi- 
tional perturbation theory. We try to explain the most 
important points of the calculation in detail. The com- 
prehensive review of the diagrammatic technique for dis- 
ordered systems can be found in Ref. [l[ We start by a 
brief discussion of the case of a weak, short-range inter- 
action potential. Although this case is artificial and is 
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unrelated to any experiment, it is governed by the same 
physics as the general problem, and it is simple enough 
to allow a transparent presentation. To generalize to 
stronger coupling, we need to recall the basic ideas of 
the Landau Fermi liquid theory and to identify the soft 
modes in the s ystem . Then we present the calculation 
leading to Eq. (2.16). Finally, to establish the relation of 
our results to existing literature, we briefly discuss scat- 
tering on a single impurity (this discussion is completely 
analogous to the one in Section || but uses the language 
of diagrams). 



A. Hartree-Fock considerations. 

The static conductivity of a system of electrons is given 
by the Kubo formula 



lim Re 



(3.1) 



l/T 



— / dT{TrMr)jf3me 



where ja{T) is the operator of the electric current at 
imaginary time r and the analytic continuation of the 
function defined at Matzubara frequencies = 27rTn 
to function analytic at Imio > is performed. 



o 
o- 




<^,. '^(^,. -a>' 



FIG. 9. Interaction correction to conductivity in the low- 
est order of perturbation theory. Here solid lines correspond 
to Matsubara Green's functions —G{ien\fi,r2) and the wavy 
line represents the interaction potential, ~V(ri ~ r2). 

Assuming that electrons interact by means of a 
weak, short-range interaction (range shorter than 
vprnin^T^ l/T), V{r) it is sufficient to consider the lowest 
order of the perturbation theory. The perturbation the- 
ory can be conveniently expressed in terms of Feinman 
diagrams. The lowest order diagrams for the interaction 
correction to the conductivity are shown on Fig. The 



Hartree term corresponds to the diagrams "a" , while the 
Fock contribution corresponds to diagrams "b" . Eval- 
uation of the correction consists of two main steps: (i) 
analytic continuation to real time, and (ii) disorder av- 
eraging. While these two steps can be performed in any 
order without affecting the result, it is more convenient 
(for technical reasons) to start with step (i) 



Although analytic continuation in Eq. (3.1) is now a 
textbook task, we include a brief discussion of the stan- 
dard procedure in the Appendix to make the paper self- 
contained. After the continuation any physical quantity 
is expressed in terms of exact (i.e. not averaged over dis- 
order) retarded and advanced Green's functions of the 
electronic system, which are defined as 



(3.2) 



where j labels the exact eigenstates of the system and 
are the exact eigenvalues, counted from the Fermi energy: 



2m 



Here U{'r) is the disorder potential. 

The resulting expression for the correction to the sym- 
metric part of the conductivity (the HalLconductivity will 
bo-discussed in a separate publicationEll) can be written 



Sa, 



a 13 



xlm' 



S^2 



-—- i I coth — - 
dn V 2T 



d^r^d^TA (3.3) 



+{«--/3}j|, 

where the extra factor of 2 in the Hartree term is due 
to the summation over electron spin. Here we denoted 
products of four Green's functions as Bp(^H)- For the 
Fock term we have 



(17;r3,r4) 



d^Tid^rc^ 
V 



;{ jrGf5(e)jfG^3(6)Gf4(e - ^U^) (3.4a) 

+J'IGU^)jIgU^)GU^ - f])Gfi(e) (3.4b) 

+2Jf Gf3(e)G||(e - n)P^GU^ - ^U^) (3.4c) 

-jrGf5(e)jf G^3(6)Gf4(e - O)Gfi(e) (3.4d) 

- Jf Gf3(e)G||(e - n)J^,GU^ - f^)Gfi(e)}, (3.4e) 
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where V is the area of the system. Equations ( 3.4c ) and 
( |3.4e| ) c ome from the diagram "b3" on Fig. ^ and the rest 
of Eq. (3.4) correspond to diagrams "bl" and "b2". For 
the Hartree term the expression is similar, 



x{j^G^,ie)J^Gt,ie)GUe-n)Gi{e) 



+J^GUe)J^Gg{e)Gilie - n)G^,ie) 



+2Jf Gf3(e)G« (e - n)J^Gi,{e - n)Gi,{e) (3.5c 



J^GUe)J^Gi,{e)Gi\(e - n)Gt,{e) 



(3.5a) 
(3.5b) 



(3.5d) 



Gf3(e)G« (e - n)J^G^,ie - f^)G^i(e)}. (3.5e) 



diagrams for averaged quantities can be constructed us- 
ing the four "building blocks" (we use the momentum 
representation since translational invariance is restored 
after averaging): 

(1) the average electronic Green's function (denoted as a 
solid line; there should be no confusion with the previous 
use of the solid line for exact Green's functions before 
averaging), which in momentum space can be written as 



(G«(^))(fc,e) 



1 



(3.8) 



2r 



(2) the disorder potential, which is assumed to be Gaus- 
sian with the correlator 



{U{n)U{r2)) 



1 



2tti/t 



In the diagrams this correlator is represented by the dot- 
ted line; 



Again, Eqs. (3.5c) and (3.5e) correspond to the diagram 
"a3" in Fig. |9|. The current operator is defined as 



eA{r} 



(v/i) /2 - (/1V/2) 



/l(^1/2(r). 



(3.6) 



In the above expressions terms corresponding to dia- 
grams "b3" and "a3" on Fig. |^ allow for at least one 
of the spatial integrations to be performed with the help 
of the identity 



dr,Gg{e)J^Gl^,{e) = -ie{n - r4fG§,{e). (3.7) 



Now it is clear that Hartree terms Eqs. ( ^.5c ) and ( 3.5c ) 
vanish identically, since there the identity ( |3.7[ ) should be 
applied with coordinates fs and being equal to each 
other. In the Fock terms Eqs. (3.4e) and ( |3.4c ) one needs 
to further multiply the result of Eq. (3/7) by the interac- 
tion potential V{r3 — r4). In the case of the short range 
interaction potential this also gives vanishing contribu- 
tion. Thus we conclude, that the diagram "a3" on Fig. ^ 
does not contribute for any form of the interaction, while 
the diagram "b3" vanishes for the short-range interac- 
tion. 

Th e sam e identity can also be applied to terms 
Eqs. (3.4cl) and (3.5d), which also vanish by the same rea- 
son. Thus the task of averaging over disorder is now sim- 
plifie d bec ause we on ly need to average two Fock terms 
Eqs. ( 3.4a ) and ( 3.4b ) and two Hartree terms Eqs. ( 3.5a ) 
and ( 3.5b| ) . These expressions contain only Green's func- 
tions of non-interacting electrons and can be averaged 
using the standard diagrammatic technique of the the- 
ory of disordered systems (see Ref. |] for review). The 





k - q 



FIG. 10. Dressed interaction vertex. 



k - q 



(3) the dressed interaction vertex F (q and are momen- 
tum and frequency of the interaction propagator) , which 
represents a geometric series in disorder potential shown 
on Fig. Hol; 



l/r 



S 



(3.9a) 



where we denote 



S 




(3.9b) 



(4) the averaged product of a retarded and an advanced 
Green's functions (sometimes referred to as the diffuson), 
where we have summed up a geometric series shown on 
Fig.|l| 



r 

A 

FIG. 11. Diffuson - geometric series of impurity lines 
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Using these building blocks we can average the prod- 
ucts of Green's functions as shown on Fig. ^ It is con- 
venient to write the averaged hi the momentum 
representation. The product Bp, which appears in the 
Fock term, can be viewed as a function of coordinates of 
the two interaction vertices and can be transformed to 
the momentum space as 



^'l(^^-^'^){BF{q,n)). (3.10) 



Using the expHcit expressions Eqs. (3.8) and ( |3.9D we can 
write the analytic form of the averaged Bp 



In the absence of magnetic field, B^'^^^ = 5°'^B 



}F(H), 



whic h is why we did not include the Greek indices in 
Eq. (^Hl) . 



;« A X A .„ AC- ^ A ;a A Cr^-rS A 




FIG. 12. Averaged product of four Green's functions. The 
wavy lines indicate S{pin ~ Pout + q) for the Fock contribution 
Bpiq) and 2nS{pin — Pout + q)S[npin] for the Hartree contri- 
bution Bh (g, ni , 712 ) . 

The Hartree contribution is considered analogously. 
One can write 



{BHin;n,r2)) 



d9i (192 (Pq 



(27r) (27r) (27r)2 

xBHi^;ni,fi2,q), (3.11) 



where Hi — (cos 9 i, sin 9 i) indictates the direction of the 
momentum. Then, disorder averaging of Bh{^] fii,n2,q) 
is performed with the help of the same diagrams (see 
Fig. |l^) but the expression for the vertices changed as 
indicated in the figure caption. 

Accordingly, the expression for the dressed vertex 
(3.9a), see also Fig. O, is changed to 



^H{n, nk\ q; ^) = 2TrS{nnk) 



1 



S 



Sh 5r — 1 ' 
Snin, q; fl) = iil — ivpqn + 1/t, 



(3.12) 



where ftk corresponds to the direction of the momentum 
k on Fig. [lo| . The final expression for Bh is similar to 

Eq. iKm 



-BH{fl;ni,n2,q) 



-27r(5 



d93 
S I 2tt 



r//(ni,n3)rH(n2,n3) 



(3.13) 



+ 



y 277 y 277 ^ 

+ {nin2) 



rif(ni,n3)rg(n2,n4) 
S'_f/(n3)5'ff(n4) 



SH{ni)SH{n2) 



and we suppressed the arguments q, $7 in the right-hand- 
side of the equation. 



We n ote in passi ng, th at by construction of Eqs. (3.4a) 
(F4b|) , (FSal) and (Hbj) that 



BF{n-r,r) = BH{n;r,r), 



and, therefore, according to Eqs. (3.10) and (3.11) the 
relation 



BF{n-q) 



d9i d92 



BH{^\ni,n2,q) 



must hol d [th is can be e asily verified using explicit ex- 
presions ( 3.11 ) and ( 3.13| )]. 

We are now prepared to calculate the temperature de- 
pendence of th e cond uctiv ity from Eq. (3.3 ). We substi- 
tute Eqs. ( |3.1lD and ( |3.13| ) into Eq. (Q. As we wiU see, 
the main contribution to the temperature dependence 
is provided by wave-vectors qr — max{T, [T /tY/"^) jvp- 
On the other hand the potential V{r) has a range much 
shorter than l/qx- This enables us to use the following 
approximations 



V 



V[rz - f4)e*«>"'^^^"'^) w 1/(0); 



d'^r^d ta 
V 



V \2kp sin 



nin2 



where V(k) in the right-hand-side of the above equations 
denotes the Fourier transform of the interaction poten- 
tial. 

Altogether, we now write the conductivity correction 

as 
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(3.14) 



d02 
2tt 



Vo 2kp sin 



nin2 



BH{^\ni,n2,q) 



Evaluating this integral (where we only k eep the temper- 
ature dependent part, see Section IIIF for details) one 
arrives to the same result Eg. (2.14), but with the co- 
efficient in the form Eq. ( ^.15 ), in agreement with the 
discussion of Section ||. 

Let us now turn to the case of the Coulomb potential, 
where the scheme of the calculation (as described so far) 
breaks down. In the Pock term we have V{Qi), which di- 
verges for the Coulomb interaction (V{q) ~ 1/g). To ob- 
tain meaningful results one needs to take into account the 
effect of dynamical screening. The Hartree term seems 
to work better since using just the static screening makes 
the result finite. However, this is wrong also, since in this 
case diagrams with extra interaction lines do not contain 
any smallness (see e.g. Fig. there the correction is 
'-^ V{2kp)V{Q)). Thus one can not justify the perturba- 
tion theory in the interaction potential. The way out of 
this problem is the standard theory of Landau Fermi liq- 
uid, which we briefly discuss in the following subsection. 





FIG. 13. (a) Si ngle i mpurity contribution to the Hartree 
term, see Section [II F for a detailed discussion (b) Second 
order correction to the Hartree term (a). 



B. Soft modes 

As we already discussed, the main contribution to 
the temperature dependence of physical quantities comes 
from the processes characterized by spatial scales much 
larger than the Fermi wave-length \p. Therefore, there 
is a scale sepacation in the problem; all the Fermi liq- 
uid parameterscil Fi are established at small distances 
of the order of A_f, and are not affected by disorder if 
the relation epT ^ 1 holds. On the other hand, all the 
temperature and disorder dependence is determined by 



infrared behavior of the system where Fi can be consid- 
ered as fixed. 

Therefore, our first step is to identify the terms in the 
interaction Hamiltonian, which may produce the biggest 
contributions at temperatures much smaller than the 
Fermi energy. This procedure contains nothing new in 
comparison with the standard identification of singlet, 
triplet and Cooper channels, see Ref. and we present 
here the main steps to make the paper self-contained. 

The original interaction Hamiltonian has the form 



Hi, 



V{q) 



'>PL (pi)V^L (P2)'0<T2 (P2 + q)i^ai {pi~q), 



and we imply summation over repeated spin indices. Soft 
modes of the system correspond to the situation when 
two of the fermionic operators have momenta close to 
each other. The difference of the momenta q* defines the 
scale 1/q* ^ Xp, which is the smallest length scale al- 
lowed in the theory. Therefore, we explicitly separate the 
Hamiltonian into a part that contains all the soft modes 
(first three terms) and a correction 6H, which does not 
contain such pairs of fermionic operators: 



Hi, 



Hn + H„ 



Hr, 



SH. 



(3.15) 



The explicit expressions for the entries of the Hamil- 
tonian ( |3.15 ) are the following. The interaction in the 
singlet channel (charge dynamics) is described by 



H. 



V{q) 



FP (Thri2) 



(3.16) 



\q\<q*,Pt 

X bPU (Pi) V'm {pi ~ q)] [i^t^ {p2) ipa2 {P2 + q)] , 

where Ui = Pi/\pi\, the dimensionless parameter 
FP (jiin2^ is related to the original interaction poten- 



tial V{q) by 



FPi0) = --V{2kpsin-], 



(3.17) 



and v is the thermodynamic density of states of non- 
interacting electrons (introduced here to make F^ dimen- 
sionless). 

Interaction in the triplet channel (spin density dynam- 
ics) is governed by 



H„ = 



\ E E 



F' 



(3.18) 



Pi 

\q\<q* 



[V'ii (Pl) Cri^^J^a^ (Pl - q)] [V'ia (P2) ^i^a^tpa^ iP2 + q)] 

where parameters F'' (jiin2^ 



are 
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F" (9) - 
Finally, the Hamiltonian 



-y 2fc.sin- 



(3.19) 




("1^2) 



X [^-.3 (P2) <.,^.4 (9 - P2)] (3.20) 

- E \ ^ [<(pi)<,,<(g-pi)] 



V-aa (P2) V'<T4 (9 - P2) 



describes singlet, F'^, and triplet, F°, pairing fluctua- 
tions. The parameters in this Hamiltonian are 



y ( 2fcF sin - j ± y ( 2kF cos - 



(3.21) 



where plus and minus signs correspond to even (e) and 
odd (o) pairing respectively. Here (T^io-a ^-re the elements 
of the Pauli matrices in spin space 



1 

1 



-I 

1 



1 

-1 



and — a^a-' . 

Deriving Eqs. ( 3.16 ) ~ ( [3.2l| ), we used the condition, 
q* <C kp- This condition allowed us to make the follow- 
ing approximation 

(Pi-P2)^«44sin2 ('^]. 



We also used the identity 



3 

(74(72 



'^(71(73'^ (72 (74 (71(73 V la 



So far, the representation (|3.15|) of original interac- 
tion is exact. The only advantage of this representation 
is that it explicitly separates the term 5H which docs 
not contain coupling to the low energy excitations of the 
fermionic system. Therefore, the contribution of 5H to 
physical quantities is regular and not infrared divergent 
[like [T /vpq^y']- Therefore, for the electron system with 
weak short range interaction, 5H can be disregarded at 
all. 

Moreover, even if the interaction is not weak or long 
range, 5H can be treated in all the orders of perturbation 



theory without generating a soft mode. If this term does 
not break the translational symmetry at short distances, 
its only effec t is to r enorm alize the i nteraction parameters 
F's in Eqs. ( |3.16|) , (|3.18D and (|3.20D and the Fermi veloc- 
ity in the non-interacting part of the Hamiltonian. For 
instance, one obtains for the two-dimensional electron 
gas with the Coulomb interaction V{q) = / {K\q\) 



FP {9) = F" {9) 



1 



where 



2 V2\ sin ■ 



kHvf 



(3.22) 



(3.23) 



is the conventional parameter characterizing interaction 
strength and k is the low frequ ency dielectric constant of 
the host material. Expression ( 3.22 ) is applicable only for 
<C 1, however, keeping it in denominator is legitimate 
for small angle scattering. 

For stronger interactioturs > 1, but still far from the 
Wigner crystal instabilitycj, ^ 37, exact calculation of 
the parameters F from the first principles (as well as their 
explicit expressions in terms of r^) is not possible. Never- 
theless, to study the behavior of the system at distances 
much larger than A p , one can still disregard the term 6H 
in Eq. (pTTsh. T hen p aramet ers F are no long er bo und 
by Eqs. (p7|), ( |3T9| ), and (|^) [or by Eq. ( |I2^ ) for 



the Coulomb interaction] but rather should be treated as 
starting param e ters f or the low energy theory. The form 
of Eqs. ( I3.16D , (|3T|) andJS^) is guarded by symme- 
tries of the system: Eq. ( 3.16| ) is guarded by t ranslational 
symmetry and charge conservation; Eq. ( 3. IS ) is guarded 
by translational symmetry and symmetry with respect 
to spin rotations; and Eq. ( 3.2C ) is guarded by all above 
symmetries and the electron-hole symmetry, which holds 
approximately at low energies. 

All the consideration above essentialW-repeats the ba- 
sics of the Landau Fermi-liquid theorycil. We reiterate, 
that this theory does not imply that the interaction is 
weak; the only assumption here is that no symmetry is 
broken at small distances. 



C. Disorder averaging 

To study the interaction correction to conductivity due 
to charge and triplet channel interactions introduced in 
the previous subsection, we follow the same route as in 
the case of the short-range interaction. In particular, the 
charge channel correction is a direct generalization of the 
Fock term. We start however with the discussion of dis- 
order averaging. 

The correction to conductivity Eq. ( 0.3| ) represents the 
first order perturbation theory in the original potential 
V{q), valid when the potential is weak. For stronger cou- 
pling we make use of the effective Hamiltonian Eq. ( ^.15 ) . 
Although the diagrams for conductivity look similar to 



12 



the Fock term "b" on Fig. ^, their content is now quite 
different. First, the wavy hne now represents the propa- 
gator for one of the soft modes in Eq. ( |3.15|) . Therefore 
the expression for the conductivity Eq. ( |3.3| ) should be 
rewritten as 



(3.24) 



xlm 



V 



x(sf (r!,f3,r4) + {a^/3}) I, 



where T>^ and T>^ are advanced propagators for charge 
and triple t channels \Dt is a 3 x 3 matrix as follows from 



Eq. ( 3.18| ), see also Sec. [HE and is the product of 



electronic Green's functions given by E q. (3.4), the same 
as in the Fock term. Deriving Eq. ( |3.24 ) we assumed that 
the spin rotational symmetry is preserved, i.e. no Zee- 
man splitting or the spin-orbit interaction is present. We 
also neglected the dependence of the interaction propa- 
gators on the direction of the electron momenta. Lifting 
of those two assumptions is straightforward but it will 
not be done in the present paper. To the leading order 
in l/kpl we can average the propagators independently 
of Bp (see e.g. Ref. p. Here we proceed with averaging 
Bp and the discussion of the propagators follows. 

We have already averaged the product Bp of four 
Green's functions for the case of the short-range poten- 
tial. There the three terms Eqs. ( |3.4e[ ), ( |3.4c| ), and ( |3.4d| ) 
vanished due to the particular form of the potential. Now 
we have to take these terms into account and consider the 
full set of diagrams shown on Figs. ^ and [ij. These dia- 
grams can be ev aluated in exactly the same way as those 
in Section [II A (where we considered a subset of these 
diagrams). 

As a re sult, the averaged Bp has a form similar to 
Eq. ( 3.11 ) and can again be expressed in terms of the 
dressed vertex F [see Eq. (3.9)]. We are still interested 



in the longitudinal conductivity and thus disregard the 
Hall contribution. Thus, after averaging the correction 



Eq. (3.24) takes the form 



oo 

2 2 [ 

6c7 = -e Vp-Kv J 



—— i I coth — 



d'^q 
(2^ 



2r2 + ! \ r^viq^ 
' [ S tS^ ' 



(3.25) 



where quantities F and S are defined in Eq. (3.9) 



;a A ^r— A .fY A 



R 





A 





FIG. 14. Conductivity diagrams, group I. Diagrams (a), 
(b), and (d) were evaluated for thp short range interaction in 
Sec. Ill A. In the diffusive regimeo only diagrams (a), (d), (e). 



were considered at ui,qvF <^ 1/r. 



R 

a 






R 


AO 




J 




xlm'^ 



Bp{n,q) = - 



I?^(f7, q) + TrV^in, g)] Bp{n, q) 
2T(^0 + i)F (F2-1)t2 



53 

v2 



S 



s-^ 2 ^ vW \ s 



FIG. 15. Conductivity diagrams, group II. Diagrams (a) 
and (b) w ere evaluated for the short range interaction in 
In the diffusive regimeo only diagrams (a), (d). 



[II A 



Sec. 

(e), were considered at uj,qvF <C l/r. 



It is important to emphasize that 



Bp{n,q = Q) = 0, 



(3.26) 
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[to s ee th is one should use exphcit expressions (3J) in 
Eq. ( |3.25 )] . This property is not accidental - it is guarded 
by the gauge invariance of the system: no interaction 
with zero momentum transfer can affect the value of the 
closed loop. 

To proce ed fu rther with the actual calculation of the 
correction (3.25) we need to specify the interaction prop- 
agator. It will be done in the following two subsections. 



dsc being the screening radius, vV^l/dsc) = 1) and we 
will not take it into account. 

Second, at frequencies smaller than the plasmon fre- 
qucn cy w e can neglect the unity in the denominator in 
Eq. (3.28), which corresponds to the unitary limit, i.e. 



1 



1 s-^ 



(3.29) 



D. Charge channel 

In this section we discuss the cha rge c hannel correc- 
tion, described by the Hamiltonian (3.16). Because the 
effective interaction is characterized by the momentum 
transfer much smaller than the Fermi wave vector, the 
Random Phase Approximation (RPA) , see Fig. |l^, is ap- 
plicable. 

To simplify further considerations, we approximate the 
Fermi liquid parameter by its zero angular harmonic 



FP{d) 



(3.27) 



this approximation does not affect the final result be- 
cause of the long range nature of the Coulomb potential 
V{q^Q) oo. 

Consequently, we write the charge channel propagator 
in the form 



l^V{q)+FP 



where the polarization operator is given by 



in 



s- 



(3.28a) 



(3.28b) 



using the notation ( |3.9b| ). 

The polarization operator (B.28b) differs from the more 
standard one (used for instance in Ref. |l^) since the dif- 
fusion approximation has not been made yet. Indeed, 
expanding the polarization operator in small Q and q we 
can recover the usual diffusive form. In terms of the scat- 
tering time it corresponds to the limit Tr <C 1. We do 
not do that here since we want to calculate the conduc- 
tivity for all values of Tt. 



The form of the propagator ( 3j8 ) and expression for 
the conductivity correction (3.25) suggests that there 
could be two contributions. First, the propagator 



Eq. (3.28) has a pole which corresponds to the 2D plas- 
mon. However, the plasmon dispersion relation is 



(vFqpiyiyV{qpi) = 2n(n+- 



i.e. {vpqpiY 1^ + r) I distances larger than 

the screening radius. According to the gauge invariance 
condition (3.26) this contribution is strongly suppressed 
(by a factor of the order ofmax[T, {T/ry/^jdsc/vF, with 



Thus the original coupling V{q) as well as the renormal- 
ization of the coupling by the Fermi liquid parameter 
Eq. (3.16) does not affect the resulting propagator. In 
other words the propagator becomes universal. 

It is important to emphasize that Eq. ( |3.29 ) gives the 
upper bound for the strength of the repulsive interac- 
tion. This is guaranteed by stability of the electron sys- 
tem with respect to the Wigner crystallization, i.e. by 
the condition vV{q) + Fq > at q < q*. Therefore, we 
always have 



^0 



1 

i^ + {vV + Fp)ii n' 



so that Eq. ( 3.2S ) is indeed the upper bound for the prop- 
agator Eq. ( |3.28a ). Note, that the above condition is 
satisfied regardless of the sign of Fq. In particular, it is 
possible to have Fq < — 1 so that the so-called compress- 
ibility of the system i//(l + Fq) is negative. This fact, 
however, has nothing to do with stability of tbe, Fermi 
liquid and does not affect transport phenomenacj. 

Using the propagator Eq. ( p. 29 ) in the expression for 
the correction Eq. ( 3.25| ) we obtain after momentum in- 
tegration 



Sac = 



r dn 


d 








f^coth — 1 


j 2^ 




V 2r; 







- arctan Or -|- — + —H{Qt) In 2 

TT TTiiT 27r 



An 



l + H{nT) 
1 



arctan ■ 



1 



In 1 



(3.30) 



where the dimensionless function H{x) is defined as 
H{x) ^ 



4 + a;2 ■ 



In the frequency integral Eq. (3.30) we single out the 
first two terms as being dominant in the ballistic and dif- 
fusive limits respectively with the rest being the crossover 
function. The diffusive limit is given by 
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(3.31) 



In the opposite limit we can replace arctanfir by n/2. 
Then the integral is divergent in the ultra-violet, but that 
large constant can be incorporated in the definition of r. 
This is done as follows: 



d f VI 
dTt— — 17 coth — 



-2T + Ef coth 



Ef 
2T'' 



(3.32) 



where Ep is put for the upper limit of the integral. This 
is consistent with the approximations in momentum inte- 
gration, where one typically relies on fast convergence in 
order to set the integration limit (otherwise determined 
by the Fermi energy) to infinity and to set all momenta 
in the numerator to the Fermi momentum in magnitude. 
Since we are interested in temperatures T <^ Ep, the 
second term is essentially a temperature independent (al- 
though infinite) constant. The temperature dependent 
correction to the conductivity is determined by the first 
term. As a result 
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^ac(T.»l) = -eVy--(l.coth- 
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(3.33) 



Integrating the full expression Eq. (B.3C) we find the 
correction valid at all values of Tr, 
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where the dimensionless function f{x) is defined as a di- 
mensionless integral: 
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(3.35) 



The factor 3/8 is introduced for convenience, so that 
/(O) = 1. The integral can be evaluated analytical in th e 
two limiting cases and the result is given by Eq. (2.17). 
In the intermediate regime the integral can be evaluated 
numerically and the result is plotted on Fig. |[ 



E. Triplet channel 



In this section we discuss the correction in the triplet 
channel. Similar to the case of the charge channel, we 
need to derive the interacti on pr opagator in the triplet 
channel and then use Eq. ( |3.25| ). As follows from the 
Hamiltonian Eq. (3.18), the triplet channel propagator is 
now a 3 X 3 matrix. Apart from this minor complication, 
the propagator can be found using t he sam e RPA ap- 
proximation as the one used in Section HID. see Fig. [T^ . 




B.. 

•WA. 



FIG. 16. Interaction propagator in the (a) singlet and (b) 
triplet channel 

Similarly to the charge channel, we take the Fermi liq- 
uid coupling F'^ to be independent of electron momenta 



F'^ (9) « F^ 



(3.36) 



Unlike the case of the charge channel, this approxima- 
tion slightly affects final results (see discussion after 
Eqs. (2.16) for the drawbacks of this approximation as 
well as for its remedies). Then the matrix equation for 
the triplet propagator has the form 



per 



F 



J2 tl.k[DT]k, (3.37) 

k—x.y,z 



where i, j = x, y, z. 

In the absence of the magnetic field and spin-orbit scat- 
tering each electronic Green's function is a diagonal ma- 
trix in the spin space, and therefore 



iii{q,n)^d,kii^{q,n), 



(3.38) 



wher e Il^ {q, n) is the polarization operator given by 
Eq. (|3.28b| ). 



Altogether, using Eq. ( 3.38 ) in the equation ( 3.37 ), we 
find the triplet channel propagator as 



F^UA{n,q)- 



(3.39) 
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Before we continu e, le t us discuss the validity of the 
approximation Eq. ( 3.36 ). Consider the situation close 
to the Stoner instabiht y Fn ^ —1. In this case the pole 
of the propagator Eq. ( 3.39 ) describes a magnetic exci- 
tation in the system. In the bahistic case (g > 1//) it is 
a slow, over-damped spin wave 

-lUK{l + F^)vF\q\. 

The main contribution to the temperature dependent 
conductivity correction comes from this pole at cj ^ T. 
The corresponding typical momenta are k* ~ T/[[\ + 
Fq)vf\- Although wBpaj'e using the momentum indepen- 
dent Fq ^ it is knownEa that fluctuations in the triplet 
channel produce a non-analytic correction to the spin 
susceptibility, so up to a numerical coefficient F'^ « 
Fq (1 — \q\lkp). Such momentum dependence can only 
be neglected if fc* < fci?(l -I- F^). This translates into 
a limitation for the temperature range where the results 
listed in Section II B arc validcB: 



T<r* w {i + f^yef 



(3.40) 



At higher temperatures x > T* our theory is not appli- 
cable. 

Having discussed the validity of the approach, we pro- 
ceed with the straightforward calculation: one has to sub- 
stitute the propaga tor E q. (3.39) into the expression for 
the correction Eq. (3.25) and evaluate the integral. The 
result of the momentum integration is given by 
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where we introduce notations: 



hi{x) =H{x;l + 2F^ 



h^{x) = H{x- 1 -f 2F^) - 1 - 2F^, + (2 + iF^)H(x) 
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(3.41c) 
(3.41f) 



- - 2FSj (1 + FSY - [FSYx' + 2FS 

Here we introduce a dimensionless function H{x] y) 

1 



H{x-y) = 



{l + vY + ixyY' 
which is related to the function H{x) introduced in Sec- 



111 D| simply by H{x) = H{x] 1) 



tion 

The expression in brackets turns into its counterpart 
in the charge channel in the unitary limit {F^ — > oo). 
Its first term describes the diffusive limit described in 
Ref. [| (the formal difference in the coefficient stems from 
the difference in the definition of the coupling constant). 
The frequ ency integral is evaluated in the same way as in 
Eq. ( 3.31 ). Similar to our discussion of the charge chan- 
nel correction [s ee e.g. Eq. ( 3.33| )], we identify the sec- 
ond term in Eq. (3.41a) with the ballistic limit (which we 
discuss in more detail in the next Section). The interme- 
diate temperature regime is described by the expression 
[which appeared previously in Section II B, Eq. ( p. 16c )]: 
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where the dimensionless function t{x\ Fq) is defined as 
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Except for th e Hmiting cases [see Eq. (2.18)] the integral 
in Eq. ( 3.43 ) has to be evaluated numerically. We plot 
the result for several values of Fq in Fig. ||. 



F. Single impurity limit 

In the previous Sections we obtained the expression for 
the correction to conductivity averaged over disorder. To 
complete the calculation we needed to separately average 
the interaction pro pagat or and use the result to evaluate 
the integral in Eq. ( 3.25 ). In doing this we assumed that 
the dimensionless conductance of the system is large or 
in terms of the scattering time rEp ^ 1. We have not, 
however, assumed anything about the relative value of 
the scattering rate and temperature. In other words, the 
correction Eq. (3.25) is valid in both the diffusive Tr ^ 1 
and ballistic Tr 3> 1 limits. It also describes the cross- 
over behavior at intermediate temperatures. 

The temperature behavior of the interaction correc- 
tion in the limiting cas es ca n of course be obtained from 
the general result Eq. ( 2.16| ). As we pointed out in Sec- 
tion II B , in the diffusive limit our results coincide with 
the standard theory, Rcf. |l|. On the other hand the 
correction in the ballistic limit is subject to conflicting 
claims in litcraturcQO. Unfortunately, neither result is 
completely correct. Therefore we discuss the ballistic 
limit in some detail, starting with diagrams before av- 
eraging (i.e. diagrams on Figs. This way we are able 
to point out exactly which diagram produces the domi- 
nant result and which diagrams were missed in existing 
theories. 

We begin by discussing the Hartree term. This contri- 
bution was considered in Ref. ^ in the framework of the 
temperature dependent dielectric function. The physical 
idea was that electrons tend to screen the charged im- 
purities and thus modify the scattering rate. In what 
follows we show which diagrams describe this process 
and how to calculate the resulting correction, which ap- 
pears to be the same (up to a numerical factor miscal- 
culated in Ref. H see below for detailed explanation). 
The important difference between the two approaches is 
that the impurity screening picture described only the di- 
rect (Hartree) interaction, while missing on the exchange 
part. The latter was later considered in Ref. p7|. We think 
that this consideration is erroneous, and we discuss the 



Fock term in Section HI F 2 



1 . Single impurity limit for Hartree term 

The goal of this discussion is to show which diagrams 
correspond to the ballistic limit of the Hartree term (as 
discussed in Section ||) and how it relates to other inter- 
action corrections we discuss in this paper. 

The Hartree term corresponds to averaging the two di- 
agrams on Fig. P, where the wavy line represents a weak 



interaction potential. In this case the diagram "a3" of 
Fig. ^ is equal to zero even before the averaging (as a to- 
tal derivative) and we only need to average the diagrams 
"al" and "a2". The rigorous procedure would involve 
dressing the interaction vertices according to Fig. and 
adding diffusons Fig. |ll| as it was done in Section~ |lI A 



(see Fig. 12), evaluating the resulting expression, and 
finally taking the limit Tt —^ oo. However, the same re- 
sult can be obtained by making the expansion by notic- 
ing that impurity line brings smallness l/Tr. Therefore, 
high temperature limit may be studied by considering di- 
agrams on Fig. |l^ directly. Such approach is completely 
equivalent to that of Ref. |6[ The result [whic h can also 
be obtained from the general expression Eq. (2.16c)] is 
similar to the one obtained in Ref. |^ (the difference is 
the extra factor of In 2 found in Ref. due to an error 
in this reference, which consists in putting the energy 
of the scattered electron on the Fermi shell rather than 
integrating over it): 



San 



-4cr 
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(3.44) 



(for weak coupling 'D{2kp) = —V{2kp) The factor of 4 in 
Eq. (3.44) can be interpreted as a result of a summation 



over four spin configurations. Although correct for weak 
coupling, this factor should be modified when stronger 
interaction is considered, see discussion above. 
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FIG. 17. Single impurity diagrams for Hartree channel. 



2. Fock contribution 

In the similar manner one can discuss the single im- 
purity contribution to the Fock term. Again, for weak 
interaction we could simply expand the result of disor- 
der averaging for the Fock term Eq. ( 3.11 ) to the leading 
order in 1 /Tt. For Coulomb inter action we would ex- 
pand Eq. ( 3.25 ), since in Eq. ( 3.11 ) the special form of 
the delta- function potential was utilized to eliminate the 
diagram "b3" on Fig. |9[ Diagrammatically, such expan- 
sions equivalent to direct evaluation of diagrams with- 
out impurity lines (but with averaged electron Green's 
functions) shown in Fig. ^ and diagrams with only one 
impurity line shown on Figs. n9, GH, and Ell 
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The evaluation of the single impurity diagrams for the 
Fock term is straight-forward and is completely analo- 
gous to the Hartree term discussed in the previous sub- 
section. The result can be written as 



(3.45) 
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FIG. 19. Single impurity diagrams for Fock channel with 
the impurity line dressing one interaction vertex. 



This result contradicts (even in sign) that of Ref. |2^. 
Here we briefly discuss the reason for this contradiction. 
We notice that one has to be careful to keep track of 
gauge invariance while evaluating diagrams for the Fock 
term. Gauge invariance manifests itself in the fact that 
any interaction at zero momentum gives no contribution 
to physical quantities, which a rc ex pressed diagrammat- 
ically as close d loo ps, see Eq. (3.26). This is indeed the 
case for Eq. (mi), where we summed up a 1 the dia- 
grams. On the other hand, any individual diagram is not 
gauge invariant. In particular, each subset of diagrams 
in Figs. p"8[]2l| i s not gauge invariant. Therefore to obtain 
the result Eq. ( 3.45| ) from these diagrams one has to dis- 
regard terms which contain higher than second powers of 
the scattering rate 1 /r. 
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FIG. 20. Single impurity diagrams for Fock channel with 
the impurity line connecting a retarded and an advanced 
Green's functions across the diagram. 
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FIG. 18. Fock channel diagrams without impurity lines. 



FIG. 21. The single impurity diagram for Fock channel 
with the impurity line connecting two advanced Green's func- 
tions. 



IV. KINETIC EQUATION APPROACH 



As we already mentioned, the contributi on fr om the 
plasmon pole is small due to the condition ( 3.26 ). How- 
ever, in Ref. |27| i t was claimed otherwise. Namely, dia- 
grams in Fig. |18| were claimed to be important for the 
plasmon correction and to give a large result, while dia- 
grams in Figs, were alleged to be not important 
for the plasmon correction. This claim explicitly vi- 
olates gauge invariance and leads to incorrect conclu- 
sions. In particular, the plasmon contribution to the 
conductivity was overestimated by a factor of order of 
(vF/dsT) ~ (Ef/T). 



Our purpose in this Section is to put the treatment 
of the interaction effects in disordered systems into the 
framework of the kinetic equation. Even though at this 
point this will not produce any further physical results, 
this proves to be more convenient for practical calcula- 
tions of more sophisticated quantities, such as the Hall 
coefficient, the thermal conductivity, energy relaxation, 
etc. The kinetic equation approach is also applicable for 
the description of non-linear effects. The main technical 
advantage of the kinetic equation is that it operates with 
gauge invariant quantities from the very beginning, un- 
like the perturbation theory, where each diagram taken 
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separately is not gauge invariant (and may produce non- 
physical divergences). 

We will presen t the final form of the kinetic equation in 
subsection IV A , and show how to operate with this equa- 
tion for the conductivity calculation in subsection IV B. 



The derivation of this equation based-pn the Keldysh 
technique for no n-equilibrium systema23 is presented in 
subsection IV C. 



point-like scattering; generalization to the finite range 
is straightforward) as well as by the self-consistent field 
generated by all the other electrons: 



Stelif} 



+ -I0 {t; e, n, r) {f{t; e, r, n))r 



(4.4) 



A. Final form of the kinetic equation 

As usual in the kinetic equation approach, averages of 
observable quantities are expressed as certain integrals 
of the distribution function f(t; e, r, n). For instance the 
averaged density is 



(4.1a) 



p{t,f)^iy J de{f{t]e,r,n))n 

— 00 

and the average current is 

00 

J(t, r) = euvF / de{nf{t] e, r, n))r 



(4.1b) 



and so on. Here v is the density of states (entering into 
linear specific heat of the clean system) at the Fermi sur- 
face and is the Fermi velocity, n = (cos 6*, sin 6*) is the 
unit vector in the direction of the electron momentum 
and angular averaging is introduced as 

de 

The Boltzmann-like equation for the distribution func- 
tion has the form 



-^^ d ^ f ^ d 
dt + vpnV + evpiftE)— — \- ujcA n ^ 

oe \ on 



/ = St{/}, 
(4.2) 



where E denotes the external electric field and lOc is a vec- 
tor with the magnitude equal to the cyclotron frequency 
corresponding to an external magnetic field perpendicu- 
lar to the pla ne a nd th e di rection along the field. 

Equations (4.1) and (4.2) neglect energy dependence of 
the velocity of electrons, which makes it inapplicable for 
quantities associated with electron-hole asymmetry, such 
as the thermopower. On the other hand, any component 
of the thermal and electrical conductivities is still within 
our description. 

All of the interaction effects are taken into account in 
the collision integral 



St{/}-Stea/} + St„{/} 



(4.3) 



The elastic part of the collision integral describes 
scattering of electrons by static impurities (we assume 



T 

The effect of the self-consistent field is described by the 
last two terms, where we introduce notations: 

/o(i;e,n,f) = - / ^IticXq ''(u;)(n,3/(t; e - w, f, n))„ 



+ 



^p + eEp — 



(/(t;e-a;,f,n))„} (4.5a) 



lf{t;e,r)^-J ^Kf{cu){f{t;e-^,r,n)U. (4.5b) 



The collision integral Eq. (4.4) preserves the num- 
ber of particles on a given energy shell: integrating 
Stei{f(t]e,f,n)} over directions of n gives zero for any 
value of e [see also Eq.(4.58)]. 



The term Ii expresses enhanced momentum relaxation 
due to static disorder. The physics of this term was dis- 
cussed in detail in Section |l|. The term /q describes 
electron scattering by non-equilibrium non-local Fock like 
potential created by all other electrons. This process is 
responsible for generation of the finite drift velocity of 
electrons. One can easily see that Iq vanishes in the equi- 
librium situation /(e, n^r) = f[e + eip(f)], V^V' = —Ea- 

The kernels Kq, Ki, and Lq entering into Eqs. ( |4.5| ) can 
be expressed in terms of interaction propagators and the 
propagators describing semi-classical dynamics of non- 
interacting electrons. Explicitly: 



Kf{uj)=lm 



(4.6a) 
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X < {D) — {n^D) - {Dn^) — {D) - (Dn^ — D) 
oqp oqp dqp 
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Here, V^{lj) — T>^{~uj) d enote s the retarded interac- 
tion propagator [see i.e. Eq. ( 3.29| )] and we introduce the 
short hand notation for the angular averaging 



{aDb} = J j^^a{n)D{n,n]UJ,q)b{n), 

3 a{n)D{n, n')b{n')D{n' , n")c{n") 



(aDbDd 



(2^)2 

dodo'de" 

{2-Kf 



for arbitrary functions a,b. The function D{n,n' ]uj,q) 
describes the classical motion of a particle on the energy 
shell tp in a magnetic field: 



iu + ivpnq + [ n x 



d_ 

dn 



D{n,n'-u:,q) (4.7) 



+ —iD{n,n']UJ,q) — {D{n,n';uj,q))nj — 27r(5(nn') 



As we have already mentioned, the elastic part of the 
collision integral is nulled by a distribution function of 
the form f[e + eip{f)] for an arbitrary /. It is the in- 
elastic term that is responsible for establishing the local 
thermal equilibrium and it has the standard form 



St^„{/} ^ J dioj deiA{Lo)f{ei) [1 - /(ei - to)] (4.8a) 
x{-/(e) [l-/(e + c.)] + [l-/(e)]/(e-c.)} 

/(e) = (/(t;e,f,n))„. 

The kernel A{u!) describes matrix elements for inelastic 
processes in both ballistic and diffusive limits. The ex- 
plicit expression for this kernel is 



Aico) 
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[Re{D)]'\D"iu,q)\\ (4.8b) 



where {D) is given by the solution to Eq. (4.7) averaged 
over angles. 

The above equations are written for the interaction in 
the singlet channel only. In a situation where both triplet 
and singlet channels are present, but the distribution 
function does not have a spin structure (no Zeeman split- 
ting or non-equilibrium spin occupation present), one has 
to replace 



V' 



TrV'r, 



in Eqs. (4.6) and 
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(4.9) 



(4.10) 



Equations (^^) - (^^) constitute the complete system 
of transport equations with the leading interaction cor- 
rections taken into account. They may be used to study 
both linear and non-linear response. We reiterate that 
they do not include effects of electron-hole asymmetry 
and in this form can not produce finite thermopower. 
The Hall effect, the thermal conductivity, and energy re- 
laxation, however, areiacluded and will be studied in a 
subsequent publicationllj. In the following Subsection we 
apply the kinetic equation approach to study the inter- 
action correction to the conductivity at intermediate and 
low temperatures and reproduce the results obtained in 
Section III by means of diagrammatic technique. The 



reason for doing so is to show how the kinetic equa- 
tion works and to demonstrate explicitly that both ap- 
proaches are equivalent. 

Closing our description of the structure of the kinetic 
equation, we discuss the range of its applicability. Any 
kinetic equation implies that the distribution function 
changes slowly on the spatial scale of the Fermi wave 
length \p and on the time scale l/tp- In our case, the 
conditions are more restrictive. First, in the interaction 
correction to the elastic collision integral we take into ac- 
count only the effect of the interaction on the zeroth and 
first angular harmonics of the distribution function. This 
implies that the equation gives the correct description for 
the interaction effects on the conductivity and diffusion, 
whereas it is not correct for description of the quantities 
involving higher angular harmonics. Second, we made a 
gradient expansion in Eq. (4.5a) and only took into ac- 
count terms linear in the electric field. This implies that 
the distribution function changes slowly on the spatial 
scale Lt = min^hvp /T, vpifi-T /TY/"^)^ and on the time 
scale of the order of Ti/T. The electric field expansion 
is justified by the condition eELx T. One can check 
that both these conditions are satisfied, if the energy re- 
laxation time is much longer than the time for the elastic 
collisions. We also did not include quantum effects of the 
magnetic field. This is justified at uj^ ^ max{T/h,T^^). 

Finally, the interaction part of Eq. (4.4) is calculated 
in the first loop approximation. It means, that it has to 
be considered as the first order correction to 1 /t. If one is 
interested in the next order interaction correction to the 
clastic part, one should take into account the second loop 
correction, which is not considered i n th e present paper. 
On the contrary, the inelastic part (4^) can be consid- 



in Eq. (4.8b) 



ered in all orders to find the zero angular momentum part 
of the distribution function; the only assumption here is 
the validity of the Fermi liquid description at energies 
smaller than ep. 



B. Conductivity calculation 

In order to calculate the conductivity at zero magnetic 
field = we look for the solution of Eqs. (4.2) - (4.4) 
in a form 
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f{n,e) = fF{e)+ftrie), 



(4.11) 



where fpi^) — l/(e^/-^ + l) is the Fermi distribution func- 
tion (aU the energies are counted from the Fermi level), 
and r is the quantity to be found and it is proportional 
to the electric field. 



We substitute Eq. ( [l.tlj ) into Eqs. (^), ( |4.5| ) 

and (4.8). The inelastic part of the collision integral [see 
Eq. ( |4.8| )] obviously vanishes, as effects of the heating are 
proportional to at least the second power of the electric 
field. As a result, we obtain an equation for F: 
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We solve Eq. (4.12) by iterations. As usual for kinetic 
equations, the solution is expressed in terms of the un- 
perturbed distribution function /j? (e) and the kernels, 
which in this case are given by Eq. (4.6): 



ra(e) 



-evpTEa 



de 



(4.13) 



+4:evFT I — 
Zn 



Kf{u)fFie~co) 
+Kf{u:)fF{e) 



de 

dfF{e - u})' 



de 



-ifF{e) I ^Lf{cj)eE,:l-fF{e^. 



Substitut ing E qs. ( |4.13| ) into Eq. ( ^.llD and the result 
into Eq. ( ^.lb| ), we integrate over e and find the conduc- 
tivity 



Sa_ 

0"D 



(T = ajj + Sa, 

dio d f , w \ 
— -TT- ^ coth — - 
Ti duj\ 2TJ 



(4.14a) 



Lpjuj) 
vft 



(4.14b) 



where the Drude conductivity is an — e^h'VpT/2. Here 
we used the fact that in the absence of the magnetic field 
all the kernels are diagonal, K"^ = Sa/iKi, Lq^ = SapLo. 
We also used the identities 



oo 

2 J defF{e) 



de 



d f , \ 
—- uj coth — - 
doj V 2T) 



dujKiiuj) = / duiLojuj) = 0. 



In order to derive explicit expressions for the kernels 
Ki and Lq we have to solve Eq. (4.7) for the function D 
in the absence of the magnetic field. The result can be 
written as 



D{n,n';uj,q}— 2TrS{nn')DQ{n,LU, if) 



-DQ{n,u,q)DQ{n',u,q) 



(4.15) 



C 



Ct- 1' 



where Dq denotes the solution of Eq. (4.7) without the 
angular averaged term (and in the absence of the mag- 
netic field) 

Do{n,uj,cf) = — — . ^ ^ ^ , ■ 
Here we used the short-hand notation 



which is similar to the n otati on 5* used in Section III [in 
fact, C = 5^ see E£^^(|j^)]. Substituting Eq. ( |4.15| ) 
into Eqs. ( 1.6a| ) - ( 4.6c ) and performing the angular in- 
tegration we arrive to 



i^i(^) = -Im|^P«(c.,q) 
'C -{-iuj + l/r) 



C-1/t 
qdq 



C - {-iuj + l/r) 
C{C-\ItY 



(4.16a) 

C - i-iLU + 1/t) ] 
C(C-1/t)2 j' 

(4.16b) 



(C-i-iu + l/T))^ 1 



C(C- 1/r) 



^o(^) 

VfT 



-Im/$P«(.,,) 



(4.16c) 



Vpq^ 



9 9 

vpq 



2rC3(C-l/r)2 [C - I/t)^ 



Together with the cond uctivi ty correction Eq. ( 4.14b) the 
above expressions Eq. ( 4.16| ) are identical to Eq. ( 3.25| ) 
obtained in Section IH by means of the standard per- 
turbation theory. Thus the kinetic equation approach is 
completely equivalent to such diagrammatic calculation. 

Integration over the wave vector q requires the 
kno wledg e of the interaction propagator. Substituting 
Eq. ( 3.29| ) for the singlet channel and Eq. ( ^.3S| ) for the 
triplet channel and performing the straightforward inte- 
gration we arrive to the results in Section II B. 
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C. Derivation of the kinetic equation 



In this s ection we derive the kinetic equation discussed 
in Section IV A. For simpHcity we show the derivation 



for the case of the singlet channel interaction Eq. (3.16). 
The case of the triplet channel can be treated in the same 
manner with minor differences (introduction of extra spin 
indices) described in the end of this Section. To keep the 
discussion at the same level as in Sec tion pil| , we treat the 
Fermi liquid parameter in Eq. ( 3.16| ) as a constant, 
simil ar to our treatment of the triplet channel in Section 

inE|. 



1. Keldysh formalism 

Here we summarize the results originally obtained by 
KeldyshEa that enable us to calculate correlation func- 
tions for any non-equilibrium distribution. 

Let us first consider a Green's function of electrons be- 
fore disorder averaging. The electron- electr on interaction 
is described by the Hamiltonian Eq. ( 3.16 ). In the path- 
integral formulation it can be decoupled from fermion 
operators using an auxiliary bosonic field 0. Then the 
Green's function can be written as 

G{xi,X2)^ I [D(b]G{xi,X2\cl))e-'^^^^\ (4.17) 



with the action defined as 



SBm^ J dtdV|i0^V'f^30|+ilogZ[<^], (4.18) 



where — Vq is the (bare; following Eq. (3.16) Vq — V{q) + 
Fq /i^) electron-electron interaction propagator and Z is 
the partition function, 



Z[(t)] = (Tce^'^^'"^''''!) 



(4.19) 



In the above expressions all the fields are defined on the 
Keldysh time contour shown in Fig. In particular, 
the fermionic fields ip^ and tjj (as well as the bosonic field 
6) can be treated as doublets 



(4.21) 



where we adopt the notation that fields with a — (-I-) 
subscript (also referred to by Greek letters in this Sec- 
tion) reside on the lower (upper) part of the contour on 
Fig. ( p^ . The time dependent fermionic operators ip are 
taken in the interaction represntation 



where Hi is the one-electron Hamiltonian which includes 
the static disorder potential as well as external fields. 

Consequently, the Green's function in Eq. (4.17) is a 
2x2 matrix. Ti me ordering along the conto ur is de- 
noted in Eq. ([4.19|) by Tc. Matrices 7" in Eq. (|4.20| ) are 
defined as 



7 = 



-1 




7 = 




1 



The Green's function G{xi, X2\4>) in Eq. (4.17) is given 



by 

G(xi,X2|0) = -L-(rcV^t (^^)^ (^2)e-*^-[^'^l). (4.22) 
Z[(j)\ 



Here, as well as in Eq. ( 4.19 ) the angular brackets (. . .) 
denote quantum-mechanical averaging. In this section 
we will use the short hand notation 

Xi = {ti, r^). 



The bosonic action Eq. (4.18) can be treated in the 
saddle point approximation: 



-iF\i 



F[q^]=F[<j> = 0] + -<J>''flcj, + O{cb'), 



(4.23a) 
(4.23b) 



(4.20) 



where CTz = diag(— 1, 1) is the Pauli matrix in the Keldysh 
space. 



Fig. pa The Keldysh contour 



where H is the electronic polarization operator, defined 
as 



Ilaf3{xi,X2) 



S^F 



S4iaixi)6(t)i3{x2) 



(4.24) 



0=0 



The quadratic expansion in Eq. (4.23b) is justified, pro- 
vided that the fields are slowly changing on the scale 
much larger than Xp. 



Let us now average the Green's function Eq. ( |4.17| ) 
over disorder: 

{G{xi,X2))d,s= I [I?0](G(^i,X2|0))d..e^'<^«[^l><'-, 



(4.25) 
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where (. . .)dis hereafter denotes averaging over disorder. 

He re w e average the electronic Green's fu nctio n 
Eq. (4.22) separately from the bosonic action Eq. (4.23). 
This approximation means that we neglect correlations 
betw een m esoscopic fluctuations of the polarizabil ity in 
Eq. ( 4.23b ) and the fcrmionic operators in Eq. ( 4.22| ) 
(which describe the motion of conduction electrons). 
This is the same approximation we used in Section |lll[ It 
is justified by the well known fact that mesoscopic fluc- 
tuations are smaller than average quantities by a factor 
of the order l/{Ki^f. 

It is convenientEil to rotate the Keldysh basis as follows 



G — ^ —o'x 
2 



1 -1 
1 1 



G 



1 

-1 



(4.26) 



In the new basis the Green's function Eq. ( [4.22 ] has the 
form 



G{xi,X2\(t)) 



fG''{x,,X2\ 
[G^{xi,X2\ 



G^(X1,X2| 
G^(xi,X2l 



(4.27) 



After the averaging over t he b osonic field and ov er the 
disorder according to Eq. (4.25) the entries in Eq. (4.27) 
acquires the following meaning/ where after integrating 
over the bosonic field the diagonal elements G-^^"^^ be- 
come the retarded (advanced) Green's functions of the 
electron system 

(G^(ti,t2)) -^ry(^2-^l)(VX^l)^^(^2) + ^H^2)^(^l)), 

where r]{t) is the Heaviside step function. The lower di- 
agonal element vanishes due to the causality, 

(G^(ti,i2)) =0, 

even before the disorder averaging. Finally, the upper 
off-diagonal element (the so-called Keldysh Green's func- 
tion) is related to the one particle density matrix 

{G''{h,t2))^~im,)i^Ht2)-^PHt2mti)), (4.28) 

the quantum mechanical averaging is performed with an 
arbitrary distribution function to be found from the so- 
lution of the kinetic equation. 

The bosonic field in the rotated basis has the two com- 
ponents 

<Pl{2) = ^ (0+ ± 0-) 

which are described by the propagators 



(4.29) 



The coupling Eq. (4.20) between the fermionic and 



bosonic fields in the rotated basis has the form: 



t f <Pl <P2 
t>2 4>1 



(4.31) 



The propagators Eq. ( [4.3C| ) are solutions of the Dyson 
equations 

P(l, 2) = Po(l, 2) + / d3d4Po(l, 3)n(3, 4)P(4, 2) 



2? - I 



n 



n-^ 



(4.32) 



and wc introduced the short hand notation (i) = (ti,ri). 
The bare interaction propagators are 



= V^^ - [V{ri ~ V2) + F^diri - r2)] d{h - ^2 



(4.33) 



Any classical external field takes identical values on the 
two branches of the contour and, hence, in the rotated 
basis has only a diagonal compo nent . 

The matrix Green's fimction ( 4.27 ) satisfies the equa- 
tion 



2m 



(xi) (4.34) 



-U{fi) - ipext{xi) I G{xi,X2\(t>) = i5{xi - X2) 

where U{r) is the potential due to the static disorder, 
Aext{xi) and (pext(xi) are the vector and scalar poten- 
tial due to the external electric and magnetic fields. 



eE ^ dtAext-y^fiext: eB^ — V x vlext (4.35) 

c 



The equation (4.34) is the basis for the further considera- 
tion. One can perform the disorder average in Eq. ( 4.34 ) 
in the leading in 1/ {Epr) approximation, which amounts 
to summation over all the non-intersecting impurity lines 
one obtains 



(0i(ti,fi)0i(i2,r2)> = -P^(ti,t2;fi,r2), (4.30a) 



L(ii,ri)</)2(t2,r2)) = -V^{h,t2;n,r2), (4.30b) 



b2{ti,n)Mh,r2)) = 2^^(ii,t2;ri,r2), (4.30c) 



2{tl,n)(f>2{t2,f2)) =0. 



(4.30d) 



2m 



4){xi) 



-ipext{xi) > G(a;i,X2| 



(4.36) 



I5(xi ~ X2) + I dx'iY.{xi,xs.\(j))G{xz,X2\ 
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E(xi,a;2|0) = ^G(a;i,a;2| 

ZmyT 



The equation (4.36) allows for semi-classical treatment 
introduced in Refs. and described in great details 

m Ref. 1^. Since we have already averaged the equa- 
tion of motion over disorder, the semi-classical approxi- 
mation now amounts to averaging the Green's function 
G{xi, X2\4>) over the distance from the Fermi surface. 
This is done in two steps: 



(4.37) 



and we will suppress the coordinate and the time ar- 
guments unless otherwise is stated. A product of such 
matrices should be understood as a matrix product in 
Keldysh space and a convolution in time: 



(4.42) 



cso 



and solutions of the homogeneous equation (4.39) are 
subject to the constraints 



f=fi-f2; R=-{ri+r2) 



dt Tr g(t,t;n,r) ^ 0, (4.43) 



1 r 



A, 



){ti,t2]n,r) 



(4.38) 



- I d^G[ ti,t2;n 



Pf -\ 

VF 



Since we follow the avenue of Ref. we will skip fur- 
ther intermediate steps, and use the semiclassical equa- 
tion written in the next subsection. 



2. Eilenberger equation 

The dynamics of the electron matrix Greenisj function 
is then described by the Eilenberger equationcj: 



d 



9{9)n - {g)ng 



2t 



(4.39) 



where angular averaging is defined as before 



de 

2^ 



n = cost/, smc 



and the covariant derivatives in Eq. (4.39) are defined as 
9*5 = dt,g + dt^g + iip{ti,r)g ~ ig(p{t2,r), (4.40a) 



Vg^Vg + iA{ti,r)g ~ igA{t2, f). 



(4.40b) 



where 



5ij5{ti -t2). 



The scalar and vector potentials in Eq. (4.40) have the 
following structure in the Keldysh space 



^'^"l A^x,{t^r) 



(4.44) 



(f>{t, r) 



ipextit,r) + (t>i{t,r) 4'2{t,r) 

(t>2{t,r) (pext{t,r) + <i)i{t,r) 



where ^pext and Aext are the external (classical) poten- 
tials due to the electric field £', 



eE = dtA, 



(4.45) 



acting on the electron system, and (f>i^2{t,r) are the aux- 
iliary fluctuating fields decoupling the interaction in the 
singlet channel. Because the singlet channel describes 
processes with small momentum transfers (smaller than 
q*, see Section IIIB), the fields (j)i,2(t,r) vary slowly on 
the scale of the 1/q*. 



The condition ( 4.30d ) enforces causality of the physical 
response functions. It is worth noticing that the decou- 
pling of interaction can be performed also using a fluctu- 
ating vector potential; our choice is strictly a matter of 
taste. 

In this formalism any observable quantity described by 
one e lectr on operator 0{p,f) is given by [see Eqs. ( 4.2S| ) 
and (|l38| )] 



Oit,f) 



dO 
2^ 



0{pFn-,r) 



(4.46) 



Here 5 is a matrix in Keldysh space 
glti,t2;n,f) 



gR 

9^ 5^ ' ' 



(4.41) 



X lim 

ti— tt 



-{g {ti,t;n,r))^ + ipext{t,r} 



where (. . stands for averaging over b oth a uxiliary 
fields 01^2 fluctuating according to Eqs. (4.3C). The 
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last term in brackets is a consequence of the ultravio- 
let anomaly, and its form is enforced by the requirement 
of the gauge invariance. 

Finally, the ele ctron ic pola rizat ion operators are deter- 
mined [sec Eqs. (4.24) and (4.38)] as variational deriva- 
tives of the solutions to the Eilenberger equation ( 4.39 ): 



n^(l,2) =n^(2,l) 



(4.47) 



2^1'^^^ 



2<50i(i2,r2) 



n^(l,2) 



de 



(ii,ti;n,fi) -I- Sg^{ti,ti;n,fi))^ 



27r 



3. Derivation of the kinetic equation 

Our goal now is to obtain an equation for the 
Keldysh function averaged over the fluctuating fields, 
{g^ {ti^ti\fi,'ri)) ff,. It is this quantity that determines 



physical observables, see Eq. (4.46). We will do this us- 



ing the non-crossing approximation for bosonic propa- 
gators (i.e. the first loop approximation for the collision 
integral), see Fig. This approximation is justified pro- 
vided that the resulting dynamics for the electrons (char- 
acterized by time r^) is slow in comparison with motion 
of relevant bosonic mode, Tt^ ^ 1. 

To do so, we notice that only two components of the 
matrix g are independent, and the other two are fixed by 
the constraint (4.43). For our purposes, we choose to fix 
the diagonal components 



g 



g 



(4.48) 



where the square root should be understood in opera- 
tor sense: as a sum of its Taylor series, with all aris- 
ing products hereafter being time convolutions, similar 
to Eq. (4.42). The two remaining independent compo- 
nents of the Eilenberger equation have the explicit form 



d_ 

dfi 



(4.49a) 



-i [(/-i (t 1 , r) - 01 {t2 ,r)]g^- ih , r)g'^ + ig^(l)2 {t2 , 7=0 



+^[/(.9^>n - {g'W + g^ig'U - 

zr 



on 



d 



g- = (4.49b) 

i [4>i {h , r) - {t2 , r )] g"" -iMh, r)g^ + ig'^4>2 {t2 , r) 



and we redefine the covariant derivatives Eq. (4.40) to 
include only the external scalar and vector potentials: 

dtg = dt^g -t- dt^g -f i [(fiext{ti,r) - (pext{t2,r)] g, (4.50) 



Vg^Vg + i 



g- 



Now we are prepared to derive the collision integral. 
We notice that due to causality {gz)<p = in all orders 
of the perturbation theory. We separate slow and fast 
degrees of freedom as follows 



gK = {gK)<i, + Sgx; gz = 5gz, 



(4.51) 



where 5g is the contribution fiuctuating with the auxil- 
iary fields and we calcula te it to first order in (j). In the 
same approximation Eq. (4.48|) becomes 



g" = 5{h h) Ig^'Sg^; 



g 



-5{t^-t2) + -5g^g'^, 



(4.52) 



[expansion up to the second order in 5g^ is un necess ary 
because terms of such kind vanish due to Eq. ([4.30d)]. 



We now substitute Eqs. (4.51) and ( |4.52 ) into 
Eqs. (4.49) and obtain equations governing the behav- 
ior of the fiuctuating parts 



dt -I- vpnV + LUc { ft X 



5g'~^[Sg'-{Sg%] 
= -2i<j>2{tur)S{t,~t2), (4.53a) 



dt + vpfiV + Qc[ n x -— 
on 



Sg"" + - [Sg" - {Sg"), 



= 2i(l>2{h,r)S{h -t2)-i [0i(ii, r) - 0i(t2, r)] (g^)^ 

+^ [{g'')^{6g^{g'')^)n {{g")n)^5g' {g'')A (4-53b) 
[(5^)0<Sg^((5^)n)0 - {{g'')^5g^)n{g'')^] , 



Solutions to the equati ons (4.53) should be substituted 
into the equation ( [4.49b| ) for the smooth part. Than the 
equation for the smooth part should be avera ged o ver the 
fluctuating fields 01^2 with the help of Eq. ( 4.3C ). As a 
result 



dt + vpn^ + ujc \ n y. 



d_ 

dn 



(4.54) 



= St,„{(g^)4-^Stei{(g^)0} 
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Here we separate the colhsion integrals into two contri- 
butions. The physical meaning of such separation will be 
discussed shortly. The first, inelastic part has the struc- 
ture 



= -i{[4>i(ti,r} - 4>i{t2,r)] 5g {ti,t2;n,f))^ 



Stm {(5^)0} {ti,t2;n,r) 



(4.55) The second, elastic contribution has the form 



Stei {(3^)0} {ti,t2;n;r) = - {{g^ {ti,t2;n;r))^)n - t2; n; f))^ 



(4.56) 



d0i 



dts I ^ [{{g'^{ti,t3;ni,T^)^F^{t3,t2;ni,n;f) - ({g^ {ti,t3;n,r})^F'^{t3,t2;n,ni;r)] 



dh J [F^{ti,t3-n,ni;r){{g^{t3M\ni,r))^ - F^(ti, tg; ni, n; r)((g^(t3, ia; n, r)),/,] , 

where the first term is ju st the ordinary impurity scattering and the remaining terms characterize interaction effects. 
The kernels in Eq. (4.56|) are defined as 



F"'(ti,t2;n,ni]f) = ^ / dtsiSg^ {ti,t3;n,r) [5g^ {t3,t2]ni,r) - 5g^ {t3,t2;n,r)\)^ 



F^{ti,t2;Ti,ni;r} ^ — / dt3{[Sg'' (tuh^ni^f) ^ Sg"" ih,h;H,r)] Sg'' ih,t2;ri,r})^ (4.57) 



The equations ( [4.54D , ( [4.55D , ( [4.56D , ( [4.53D , ( |f.30D , and 
(4.32) constitute a closed system of kinetic equations. Al- 
though sufficient for description of interaction effects in 
disorder systems, these equations are inconvenient for an- 
alytical calculations because the expressions for the col- 
lision integral are nonlocal in space and time. To sim- 
plify further calculations we will use the assumption that 
{g^{ti, t2\ n, r))^ is a smooth function so that a gradient 
expansion will be possible. 

Before embarking on such calculation we pause todis- 
cuss the phys ical distinction between the elastic ( 4.56| ) 
and inelastic ( |4.55|) c ollision terms. One immediately no- 
tices from Eq. (|4.5q) that 



dO Ste;(ii, ^2; n] r) = 0, 



(4.58) 



for any ti and t2. This indicates that this part of the 
collision integral preserves the number of particles on 
a given energy shell [see below for explicit connections 
betw een t ime representation and energy representation 
Eq. ( ^J5|) ]. 

The inelastic term ( 4.55| ) does not vanish after angular 
averaging. Therefore this part does promote electrons 
between energy shells. However, we notice that 



St„{(g^')0}(ii,ii;n,r) = O, 



(4.59) 



for any direction n. Taking coinciding time arguments 
is equivalent to integrating over the whole energy spec- 
trum [see Eq. ( 4.65| )], so that not only the total number 
of particles is conserved, but the total number of parti- 
cles moving along a given direction n is conserved (i.e. 
inelastic forward scattering). 



Let us now perform the ac tual ca lculation of the colli- 
sion integrals. We solve Eq. (4.53a) and obtain 



5g^{h,t2\n,r) = 2i5{h-t2) 



(4.60) 



X / dYidt3(t>2{vi,t3) I —D{t3 - ti,n ,n]ri,r) 



D{t;n,n';ri,f2) = 



72^' 



,iq{fi—r2) — i^t 



D{n, n';uj,q), 



where the diffuson propagator D is define d in Eq . (4.7). 

To simplify the analytic solution of Eq. ( 4.53b ), we as- 
sume without loss of generality that (5^)0 varies slowly 
on the spatial scale Lt = VF'min{\/T, \J t /T), and also 
a slow function of ti +t2 on the time scale ~ 1/T. These 
assumptions are consistent with the first loop approxi- 
mation we already invoked. 

In what follows we will keep only the zeroth and first 
angular harmonics (which is consistent with assumption 
about the spatial smoothness) in the direction depen- 
dence of the Keldysh function: 



(5 {ti,t2;n,r))^ w {g {ti,t2;n,r))r 



+2n{ng {ti,t2\n ,r))r 



(4.61) 



This approximation does not affect results for any rele- 
vant quantities. From now on we will suppress the ex- 
plicit sign of averaging over the fluctuating fields because 
we will not be dealing with non-averaged quantities any- 
more. 



We n ow su bstitute Eq. (4.61) into the right-hand side 
of Eq. (4.53b) and obtain 
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+Sg2{ti,t2;n,f) 



(4.62) 



The first term in Eq. (4.62) is proportional to tlie field 
and gives contributions to both the elastic and the 



inelastic parts of the collision integral. To obtain non- 
vanishing contribution to the latter we have to do each 
one of the following: (i) take into account the first angular 
harmonic; (ii) perform the first order gradient expansion; 
(iii) expand up to the first order in external fields, A^xt- 
The result is 



dt[(t)i{fi,ti-t) - (j)i{fi,t2-t)] / —D{t,fi,n';r,fi) 

X [{g^ih -t,t2-t;ni,r))n, 

+2n'{nig^{ti -t,t2- t; ni, r))si 

+ (ri - V(<?^(ii -t,t2- t; 7^1, r)),i, } 



(4.63a) 
(4.63b) 
(4.63c) 



where the covariant derivative is defined in Eq. ( 4.50| ) and we neglected higher order derivatives of the external fields. 
Expansion in the time coordinate ti + t2 (using the covariant derivative dt) is not necessary because it produces a 
negligible correction to the inelastic collision integral and does not affect the elastic one. 

The second term in the RHS of Eq. ( 4.62| ) is proportional to the field 02, and according to Eqs. (4.57) and (4.30c]) 
it does not contribute to the elastic collision integral. Therefore, it is sufficient to keep only the zeroth angular 
component and neglect gradient terms at all. This yields 



Sg2{ti,t2;n,r) 



dO' d9" 
27r ^ 



dridtD{t, n, n' ; r, fi) |2i02(fl,^i - t)S{ti - 12) 



(4.64) 



+ -(5 (^1 - t,h\ni,r))n^{g {h,t2 - t;ni,r)).fi^ [{D{ti - ts, 7I", ni; fa, fi))si - Dit^ - ig, n", 71'; fa, fi)] .(/)2(r'2, ^4) 



As we already mentioned, g{ti,t2) has a much faster 
dependence on the difference ti — t2 then on the sum 
ti + t2- Therefore it is more convenient to use a temporal 
transformation of the Green's function 



g {ti,t2;n,f) 



de 
2^-' 



tl+t2 



Je(t2-ti) 



(4.65) 

which defines the precise notion of energy e in this con- 
text. We introduce the sa me tra nsformation for the prop- 
agators of auxiliary fields ( 4.30| ) 

'ti+t2 



Juj{t2-ti) 



(4.66) 



The transformed functions have the symmetry property 
(hereafter we omit the K superscript for brevity since we 
are only dealing with the Keldysh function) 

dco , 



g{t,e) = -g{t,-e); (4.67) 
{t,Lo-ri,r2) =V^ it,-Lo;f2,ri); 
V^{t,uj;fi,f2) =V^{t,-Lu;r2,ri). 



Now, we are ready to obtain the explicit form of the 
collision integral. We start with the inelastic contribu- 
tion and perform the following three steps: (1) substitute 
Eq. ( [t.63a] ) and ( [4. 64] ) into Eq. ([4.8a| ); (2) average over 
the fields (/)i.2 with the help of Eq . ( 4.30 ); (3) perform 
the temporal transformation (4.65) of the r esult. As a 
result we obtain with the help of Eqs. (4.66) and (4.67) 
the following form of the collision integral: 



St» {g^} (i, e; f) = ~- / d\, / — (t, f, n) [{D{lu; f, n)) + {D{-uj; n,r))] 



X [(5 (t, e; n, r)),i - {g {t, e-uj;n, r))fi] 



'2^ 



dfi dr2 



duj 
2^ 



[V^ (t, Lj; n,r2) ~ {t, w; ^2,7=-!)] [{Diio; r, n)){D{-uj; f, r^,)) - {Dico; r, ri)D{-u; r, 7^,))] 



X {g {t, e + uj;n, r)).ft{g {t, e; n, r))n 



(4.68) 
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where the angular averaging of the difFusons is defined after Eqs. 

Now, we have to express th e bos onic propagator in terms of the fermionic polari zation operators. The polarization 
operators are gi ven b y Eqs. (4.47), where we now substitute Eqs. ( 4.6C| ), ( |4.63a ), and (4.64). After the temporal 
transformation (4.65) we find 



'5(n - ^2) + ;^(£'(t^;?^i,f2)) / de[{g{t,e;n,r))n- {git,e-uj;n,r))n] 



v[S{fi - ^2) + iLL!{D{uj;fi,r2)}] 



(4.69a) 
(4.69b) 



ff'(i,w;ri,f2) = ^ I dr 



{D{u;;f,n)){D{-u;;f,r2)) - {Diu;;r,n)D{-u;;f,r2)) 



(4.70) 



The last step in the calculation of the interaction propagators is to solve Eq. (4.32) with the polarization operators 
Eq. (4.69a). This gives V^'"^ in the form given by Eq. (3.28a) and for the Keldysh component we obtain 



(4.71a) 



Also we can relate the difference of the retarded and advanced propagators which enters the collision integral Eq. (4.6J ) 
to the polarization operators: 



pi?, _ ^ pi? ^jji?, _ jjAj 



(4.71b) 



To obtain the final form of the inelastic part of the collision integral Eqs. ( |4.8| ) we need to substitute Eq. (4.71) 
into Eq. (4.68), while using Eq. (4.69a) for 11''' — 11"^. In addition, we note that 



{D{uj; q)D{-LO- -q)) - {D{uj- q)) {D{-oj; -q)) = {D{uj; q)) + {D{-uj; -q)) 



Finally, we introduce the gauge invariant distribution function / as 



f{e,t;n,r) = ^ - ^g{e + Lpext{r),t;n,r) 



(4.72) 



and obtain Eqs. (4 



The calculatio n of t he elastic p art o f the collision integral is completely analogous. We substit ute E qs. ( 4.63a )- 
(4.63c) and Eq. (4.60) in to Eq s. (4.57) and average over fluctuating fields with the help of Eq. ( 4.30| ). After the 
temporal transformation (4.65) we find 



F^{e,t;ni,n2,r) = F^{e,t;ni,n2,r}* ^ ^ / ^ f dndr2V^{uj,n,r2) 

X / ^^^^[D{u]n3,n2,f2,r^ - D{uJ]n3,ni,f2,T^]D{uj;ni,n^ 



+2n4{n'g(t,e- uj;ni,r))ff 



(4.73a) 
(4.73b) 



+ {n-r-)h + ^§-){git,e~.;n 



(4.73c) 



Deriving Eq. (|j|) we use the fact that / dujV^{uj)D{uj) = 0. 
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We substitute Eqs. (4.73) and (4.61) into Eq. (4.5S). We ex pand the result into angular harmonics. The zeroth an- 
gular harmonic vanishes because of the conservation law ( 4.58 ), and the first harmonic gives (we write only interaction 
correction to the collision integral) 



duj 
2^ 



+ UaK'^'^ {u}){npg{t] e-uj,r, n))n{g{t; e, f, n))„ 



(4.74a) 
(4.74b) 



{g{t;e,r,n))n ( + ) {9{i\^ - ^,r,n))r. 



(4.74c) 



Here the kernels K and L are given by Eq. ( |4.6| ). Each labeled separately term in Eqs. (4.74) corresponds to ones in 
Eqs. ( 4.73 ) and in Eqs. ( 4.63 ) so the origin of terms can be easily traced. 



Fi nally, we use the gauge invariant distribution func- 
tion (4.72) instead of g and we arrive to Eqs. (4.5). 

Closing this section, we remark that the above treat- 
ment can be easily generalized to include other chan- 
nels as well as the higher angular harmonics of the Fermi 
liquid constant. For the latter task one has to intro- 
duce angle dependent auxiliary fields 0i.2(^, ^, i) and use 

The triplet channel requires introduction of the cou- 
pling of the form hi^2{^Tt)a, where ct^, j = x,y,z are 
the Pauli matrices in the spin space, and hi 2{r,t) are 
the auxiliary fiel ds. Accordingly, each bosonic propa- 
gator fro m Eq . ( 4.30 ) becomes a 3 x 3 matrix. The 



^0 



equation ( 4.32| ) retains the same form with the matrix 
multiplication in Keldysh and spin spaces implied. The 
equation ( [4.33| ) becomes 



f2)S{h-t2), (4.75) 



and Eq. (4.47) is modified to 



n,^(i,2) = nf,(2,i) 



^j^[5,25,.+ 



n^(i,2) = Tiv 



(4.76) 



7r(Trcr,Jg^(ti,ti;n,ri))0 
4<5M(i2,r2) 



dO (Trcri(5g^(ii,ti;r7:,ri) + Trcr,(5.g^(ti,ii;n,fi))^ 



27r 



A5hi{t2,?2) 



where trace is performed in spin space. 

Further derivation consists of a repetition of the steps 
described in this section, and in the absence of the spin 



structure of the distribution function, fi, 



suits in Eqs. (4.9) and (4.10). The spin-orbit interaction 
or Zeeman splitting by external magnetic field slightly 
changes the results, but we will postpone-the correspond- 
ing analysis until the future publicationO-3. 



Finally, the Cooper channel interaction (B.20) can be 
treated in the same manner by introducing auxiliary 



fields in the Gorkov-Nambu space. We will not discuss 
this question further in the present paper. 



V. CONCLUSIONS 

This paper is an attempt to consistently describe the 
effect of electron-electron interaction on longitudinal con- 
ductivity of disordered 2D electron gas at T ^ Ep. Our 
results are valid for an arbitrary relat ion between T and 
h/r and are summarized in Section 11 B. At low tem- 
peratures Tt <^ Ti we reproduce the logarithmically di- 
vergent Altshuler-Aronov correction. At higher temper- 
atures Tr > h, i.e. in the ballistic region, we found the 
linear temperature dependence in accord with Refs. |6|j2^. 
However, even the sign of the slope of this dependence de- 
pends on the strength of electron-electron interaction in 
cont radic tion to the results of Refs. |6|,p7| (see Sections || 
and IIIF for discussions of this discrepancy). 

We deliberately did not compare the theory with ex- 
perimental data, postponing this comparison until the 
publication of theoretical results for Hall conductivity 
and magneto-resistance in the parallel field. For com- 
parison with data obtained for Si-MOSFET samples the 
valley degeneracy should be taken into account (the de- 
genaracy may increase the numerical factor in Eq. (2.16c) 
by as much as a factor of 5 in the case of low intervalley 
scattering). We also relegate the corresponding discus- 
sion to a separate publication. 

Finally, we derived a kinetic equation to describe the 
effect of electron-electron interaction at arbitrary Tt. 
The advantage of this approach is that it turns out to 
be more convenient for practical calculations of transport 
properties in magnetic field as well as thermal transport 
properties. 
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APPENDIX: 



FI G. 23 . Integration contours for analytic continuation of 
Eqs. ( [All ). 

In this Appendix we show in some detail the procedure 
of analytical continuation that lead s to the expression for 
the interaction correction Eq. (3.3) to the conductivity in 
terms of exact Green's function of noninteracting disor- 
dered system. The structure of the current correlator 
is 



l/T 



JaGiien + inn)JpG(ien) (Ala) 
JaG{ien + iftn)^p{ien + i^n, ien)G{ien), (Alb) 



where = 7rT(2n + 1) is the fermionic Matsubara fre- 
quency, G{ien) is the exact Green's function of the inter- 
acting system (diagrams 1 and 2 on Fig. ^ are the first or- 
der correction to the Green's function), r^(je„-|-ir2„, ze„) 
is the vertex function (not to be confused with disorder 
averaged interaction vertex F from Sec. III). Diagrams 3 
on Fig. ^ are the first order correction to the vertex func- 
tion. The current operator is defined in Eq. (3.6). Note, 
that we omit the spatial coordinates and the integration 
whenever it should cause no confusion. 



We perfo rm t he analytic continuat ion i n each term 
( Ala ) and ( Alb ) separately. In Eq. ( Ala ) we use the 
standard procedure 



detanh— (. . .), 



(A2) 



We 



where integration contour is shown on Fig. t 
deform this contour to form Ci, use the facts that 
ta.nh{e + in„)/2T ^ tanhe/2r, and G{e±iO) = G'"(^)(e) 
and wc obtain 

Mi{in„) = J„G(ie„ + iVln)Ji3G{i€n) 

= J Atanh^{j„G(6 + *a„)J^[G«(6)-G-^(6)] 

+J„ [G«(e) - G^{e)] JfiG{e - ^a^)} (A3) 

In the form ( [A3| ) frequency is present only in func- 
tions which may have singularities only on the real axis, 
so that the required analytic continuation can be easily 
performed: 

Mi{Ld) = Mi(ifi„ uj + iQ) 



de 



■ tanh ■ 



|j„G^(e + L.)J^ [G«(e)_G^(e)] 



-fj, [G«(e) - G^(e)] J^G^ie-cu)} 



(A4) 



Thus o ne obtains for the quantity entering into conduc- 
tivity (O) 



lim Im 



Miiuj) 



Re J ^tanh^J^G^ie)Jpd,G^{e) 
Jo.G'^{e)JpG^{e) 



de f d e \ 



(A5) 



Equation (A5) can be further simplified for the calcula- 
tion of the symmetric part of the conductivity 



Re 



de f d e \ 

— — tanh — 



(A6) 



Jo.G^{e)JpG'^{e) + Jc.G'^{e)JpG^{e) + (a ^ /3) 



Term ( Alb ) is considered analogously. We find similarly to Eq. (A3) 
M2{iVln) = T'^JaG{ien + iVln)Vi3{ien + i^n, «en)G(ie„) 

= J Atanh^{j„G(e + zaO [r^(e + e + iO)G^(e) - r;3(e + e - iO)G^(e)] 
+Jc, [G^(e)r/3(e + tO,e- iiln) - G^{e)Tp{e -iO,e- iQn)] G{e - 01,,)} (1.7) 
Using analytic properties of the Green's function and the vertex function, we perform the analytic continuation and 
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obtain 

A/2(w) = Ahiiflr, ^ UJ + iO) 

= lim lim / — ^tanh-^jjaG-^fe + cj) rr«(e + w + i(5i.e + i(52)G^(e)-r/5(^ 

+Ja [G^{e)rfi{e + iO,e - Lo - iO) - G^{e)rp{e - i62,e - to - iSi)] G^{e-Lu)y (1.8) 
In the appropriate frequency limit, we find 



N2 = — lim Im 



lim lim V R(e ~ 182,^1 — i5i)G (ei) 

■51^0+52^0+ 



e + iO,e-iO)G^(e) 



(1.9) 



Further calculation requires specification of the form 
of the self-energy and the vertex function. We have to 
find both in the first order in interaction propagator, and 
expand Eq. (|A6D up to the first order: 



Re / — ( — tanh — 
47r V 2T 



(1.10) 



X J„G^(e)S^(e)G^(e)J^G^(e) 
+Jc.G«(e) j0G^(e)S-^(e)G^(e) + (a ^ /3) 



For brevity we consider only the "Fock" contribution of 
Fig. I b: 

S(ien)i2 = r^I?i2(if2m)Gi2(ie„ - 10.^) (l-H) 



where T) is the bosonic propagator, and we restored the 
notation for spatial coordinates. In the same order 



[r(ie„,ie 



(1.12) 



After analytic continuation similar to that in the 
derivation of Eq. (A4) we find 



S^2(e) = - / ^coth^[lmP5^2(f^)]Gf2(f-f^) (1.13) 



+i 



— tanh -n^/oA 

47r 2T 



and for the vertex function we have two cases: 



lim lim r/3(e- i(52,ei -i^i)= - / — coth — [lmX»d(ll)l [G'^(e- f7)J^G^(ei - m 
5i^o+52^o+ 7 27r 2T L JL 



dn , e — 57 „ 4 

-I / — tanh I?r9 

47r 2r ^2 



{n) {G^{e-n)~G^{e-n)^Ji3G^{ei~n) 

+i I '^ian\i^-^^Vt2{^)\G^[e-n)JfJG^{ei-n)-G'^{ei-n) 
47r 2r 
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(1.14a) 



r/3(e + «0,e-iO)= 



+i 



f dfl 


n - 


/ 7:— coth 




/ 27r 


2T . 


df] , e 


~ n 


— tanh - 


47r 


2T 


dn , e 


- n 


— tanh - 






2T 



Iml?^^2(f^) G^{e-n)Jf}G'^{e~n) 



f2(r!) [(G^(e -n)^ G^{e - 57)) JpG^{e - n) 
^^in) [G^(e - n)Jfj (G^(e - 17) - G^(e - 57)) 



12 



(1.14b) 



We now substitute Eq. ( 1.13 ) into Eq. ( 1.10 ). We use the fact that the combination containing only retarded or 
only advanced Green's functions vanish upon the disorder averaging. Moreover, the average of the combinations like 
G(e)G(e — 57i) . . . G(e — 57 ^r) does not depend on the energy e, which enable us to perform the integration over e using 

,e-57d f O 

rfetanh tanh — = —2 — 5 2 coth — 

2T de 2T dn \ 2T 



We find using 1)^(57) = [D'^{^n)] 
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sym 



Im 





■ d 








Ocoth — I 


1 8^ 


dn 





(1.15a) 



dn 

4^ 



n 



2rsinh2 ^ 



ImV^^i^) Re J„G^(e)G^^^(e - f7)G^(e) J^G^(e) + (a /?). 



(1.15b) 



The same manipulations are performed with substituition of Eqs. (1.14) into Eq. ( |l.g| ). One finds for the symmetrized 
part 



^iVf™ = Im 



dn 

8^ 



d_ 

dn 



n coth 



n 
2r 



(1.16a) 



G'^(e) J„G^(e)J \G'\t ~ n)JpG'\e - n) 



- 2 



G^(e)JaG^(e) G^(e - 17) J^G«(e - 1^) 



dl7 r r2 
3 



2T sinh" 



^ImlJ^i^Cf^)] { [G«(e) J„G^(e)] [G^(e - n)JpG\e - n)] ^ J + (a ^ /3). 



(1.16b) 



Total cor rect ion to the conductivity is just A^i + In the Hartr ee-Foc k app roximation = —V{q) a nd we 
obtain Eq. (3^V In the c ase for the stronger interaction terms ( 1.15a ) and ( 1.16a ) are added to produce Eq. ( 3.24 ) 
and terms (1.15t), (1.16t) give rise to the inelastic or so-called dephasing term: 



5a 



deph 



dn r n 

8^ 



2rsinh^^ 



Iml?^^2(^^) Re<^ 2J„G^(e)Gf2(e-17)G'^(e)J/3G^(e) 



G^(e)J„G^(e) G^(e - 17) J^G^(e - 17) \+{a^(3). 



(1.17) 



In our leading approximation in l/EpT this term van- 
ishes, see Fig. The role of this term in the temper- 
ature dependence of weak localization correction is dis- 
cussed in detail in Ref. |l^. 








FIG. 24. Cancellation of inelastic term ( 1.17 ) in the leading 
ladder approximation. 
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